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Preface 


Due  to  the  recent  interest  in  electron  beam  sustained  lasers,  the 
effects  of  the  beam  on  the  energy  distribution  of  free  electrons  in  a 
plasma  has  gained  considerable  attention.  This  thesis  attempts  to  pro- 
vide some  insight  into  the  problem.  Although  the  results  obtained  are 
a bit  brief,  it  may  provide  a starting  point  for  further  analysis. 

I would  like  to  extend  my  appreciation  to  my  advisor,  Capt  Allen 
Hunter,  whose  insight  into  the  problem  has  proved  invaluable.  I would 
also  like  to  extend  my  sincere  appreciation  to  Judy  Csizmadi  for  her 
help  in  the  preparation  of  this  thesis  and  whose  optimism  provided  the 
incentive  for  its  successful  completion.  And,  finally,  I would  like  to 
thank  Bill  Ercoline  for  the  time  he  spent  in  gathering  data  used  in  the 


various  numerical  evaluations. 
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Abstract 


An  attempt  is  made  at  determining  the  effects  of  an  external  flux 
of  high  energy  electrons  on  the  electron  energy  distribution  function. 

An  analytic  model  is  developed  for  a simple  system  and  compared  against 
numerical  results.  Normalization  problems  precluded  detailed  investi- 
gation of  the  analytic  solution. 

Numerical  solutions  were  obtained  for  a more  complex  system  using 
a Boltzmann  computer  code  developed  for  this  and  three  related  theses. 

It  was  found  that  an  electron  beam  tends  to  "mask"  the  effects  of  an 
applied  dc  field  in  the  high  energy  region  of  the  electron  distribution 
function.  This  accounts  for  an  uncharacteristic  dip  in  efficiency 
curves  for  processes  with  high  threshold  energies  as  the  magnitude  of 
the  applied  field  increases.  The  beam  also  accounts  for  an  insensitivity 
of  pumping  rate  constants  to  variations  in  applied  field  strength  for 
processes  with  high  threshold  energies. 
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THE  EFFECTS  OF  BEAM  CURRENTS  ON 
KLEc  RON  ENERGY  DISTRIBUTIONS 

T n trod  vi  c t i o n 

Gas  lasers  have  provided  plasma  physicists  with  a variety  of  inter- 
esting problems.  One  such  problem  involves  the  "pumping"  of  a laser  by 
a mono-energetic  electron  beam.  In  order  to  insure  proper  laser  design, 
it  is  essential  to  determine  the  distribution  of  energy  among  electrons 
in  the  laser  cavity  [referred  to  as  the  electron  energy  distribution 
f(e)  ] and  how  this  distribution  changes  in  time. 

The  Boltzmann  equation  is  one  of  the  most  common  tools  used  to  des- 
cribe this  process.  It  relates  the  time  rate  of  change  of  the  electron 
energy  distribution  to:  1)  energy  losses  due  to  electron  collisions  with 
heavy  particles  (ions  and  neutral  atoms),  electron  collisions  with  cavity 
walls,  electron  diffusion,  etc.;  2)  energy  gain  due  to  collisions  with 
excited  heavy  particles,  heating  mechanisms,  etc.;  and  3)  free  electrons 
being  created  via  field  driven  ionization  and  various  external  sources. 

Numerical  solutions  to  the  Boltzmann  equation  have  been  obtained 
through  Boltzmann  computer  codes.  However,  these  codes  produce  only  a 
single  solution  for  each  set  of  physical  parameters;  hence,  it  is  diffi- 
cult to  determine  the  effect  of  each  parameter  on  the  distribution  func- 
tion. Because  of  this  problem,  laser  design  is  a somewhat  trial  and 
error  process. 

Due  to  complexity  of  the  problem,  a complete  analytic  solution  has 
not  yet  been  derived.  However,  by  making  various  assumptions,  depending 
on  the  physical  parameters,  terms  in  the  Boltzmann  equation  can  be  modified 
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or  neglected  in  an  effort  to  obtain  partial  analytic  solutions.  These  par- 
tial solutions,  although  not  completely  accurate,  provide  a great  deal  of 
insight  into  parameter  effects  and  serve  as  a significant  asset  to  laser 
design. 

One  design  area  currently  of  interest  to  the  Air  Force  deals  with 
electron  beam  sustained  lasers.  This  thesis  attempts  to  provide  some  in- 
sight into  this  problem  through  numerical  solutions  and  a partial  analytic 
solution  to  the  Boltzmann  equation  with  special  emphasis  placed  on  the 
effects  of  electron  beam  current  and  an  applied  dc  electric  field  on  the 
electron  distribution  function. 

Chapters  I and  II  provide  some  background  into  the  problem  cf  tergy 
deposition  by  an  external  flux  of  high  energy  electrons,  introduce  the 
Boltzmann  equation,  and  define  most  of  the  terms  and  symbols  used  in  this 
thesis.  The  computer  code  used  to  provide  the  numerical  solutions  is  dis- 
cussed in  Chapter  III  and  a summary  of  the  derivation  of  the  analytic  solu- 
tion is  presented  in  Chapter  IV.  Chapter  V presents  the  results  of  the 
numerical  and  analytic  solutions  and  Chapter  VI  contains  concluding  remarks 
and  recommendations. 


2 


, I.  Background 

To  induce  lasing  action  in  a gas,  it  is  necessary  to  excite  a certain 

percentage  of  the  atoms  serving  as  the  laser  medium.  One  means  of  accom- 

plishing this  is  by  applying  a large  dc  field  across  the  laser  medium. 

The  free  electrons  in  the  gas  are  accelerated,  gain  energy,  collide  with 
the  laser  medium,  and  excite  the  particles.  When  the  proper  level  of  exci- 
tation has  been  reached,  lasing  action  can  occur  (of  course,  it  is  a bit 
more  complicated  than  this) . This  type  of  excitation  process  is  referred 
to  as  discharge  pumping. 

For  optimum  discharge  pumping,  a specified  number  density  of  electrons 
in  a given  energy  range  must  be  maintained.  If  too  few  electrons  are  pre- 
sent, not  enough  particles  can  be  excited;  if  too  many,  the  free  electrons 
and  corresponding  positive  ions  reabsorb  the  photons  being  emitted  when 

the  excited  particles  return  to  their  ground  state.  In  either  case,  the 

optimum  power  output  of  the  laser  will  be  severely  reduced.  Therefore,  one 
problem  in  discharge  pumping  is  to  control  the  number  of  electrons  in  this 
critical  energy  range.  This  is  a complicated  process  since  free  electrons 
are  prone  to  recombine  with  positive  ions,  and  can  lose  energy  via  elastic 
and/or  inelastic  collisions  with  heavy  particles,  other  electrons,  and 
cavity  walls.  Therefore,  the  dc  field  must  be  large  enough  not  only  to 
accelerate  electrons  to  an  energy  state  capable  of  exciting  particles,  but 
to  accelerate  them  to  a state  capable  of  ionizing  particles,  thus  creating 
a source  of  secondary  electrons  to  compensate  for  free  electrons  lost  to 
attachment  and  recombination  reactions.  These  high  energy  electrons  may 
be  relatively  ineffective  in  exciting  the  appropriate  laser  level.  While 
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t ln-y  provide  thci  nocoasa  t y .sour  oh  of  aecoudaiy  elections,  the  oveiall  **  t t i — 
elonuy  of  the  laser  may  decline. 

Kbjure  l illustrates  on  optimum  election  enemy  ill  at  t Unit  Ion  function 
foi  excitimj  the  Idsoi  mud  l um.  In  thin  case,  a hlijh  pel  cent  ai|o  of  t ho  fum 
flout  t diis  in  t ho  laser  cavity  huvo  fUfujioM  compatible  with  t ho  excitation 
cioas  sect  Ion,  u , while  low  are  available  to  piovitto  the  nocessaty 
ion l /.at  ion  to  compmrtal  o tor  election  losses. 


Kii|uie  l Optimum  election  Kneruy  P i at  i i hut  ion 


Kiijuro  - illustrates  an  elei'tion  enemy  distribution  function  in  which 
a finite  number  of  elections  are  maintained  with  enemies  compatible  with 
the  Ion  i station  cross  sect  Ion,  o , . More,  fewot  electrons  are  available 

foi  excitation  thus  teducinij  the  ovoiall  officit'ncy  of  the  laser. 
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Function 


Figure  2 F.lectron  Energy  Distribution  to 
Meet  Ionization  Requirements 

Thus,  it  can  bo  shown  that  if  the  secondary  electron  production  can 
be  accomplished  without  modifying  the  optimum  electron  energy  distribution 
for  excitation,  the  efficiency  of  the  laser  can  be  greatly  increased. 

An  alternate  means  of  producing  secondary  electrons  is  to  bombard  the 
gas  tube  with  a mono-energetic  electron  beam.  These  high  energy  electrons 
collide  with  neutral  atoms,  ionize  the  atoms,  and  produce  the  necessary 
secondary  electron  spectrum. 

Thus,  it  can  be  seen  that  the  manner  in  which  those  secondary  electrons 
modify  the  electron  energy  distribution  is  of  considerable  interest  to  agen- 
cies engaged  in  laser  research. 

To  analyze  the  effect  of  an  electron  beam  on  the  electron  distribution 
function,  a quantitative  understanding  of  the  manner  in  which  energy  is  lost 
by  the  injected  flux  of  primary  electrons  is  required. 

Electron  Beam  Energy  Deposition  (Ref  12) 

In  examining  the  energy  loss  by  the  primary  injection  of  charged  par- 
ticles, one  can  experimentally  determine  the  energy  lost  per  unit  length 
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of  penetration,  -de/dx  . Given  this  and  the  effective  ionization  poten- 
tial, I , the  number  of  ion  pairs  produced  per  unit  length  can  be  deter- 
mined. Therefore,  the  volumetric  ionization  rate,  SQ  , due  to  the  pri- 
mary beam  of  particles  is 
so  3 3eb  (-de/dx) 


q 


(1-1) 


where  j0b  is  the  electron  beam  current  density  and  q is  the  charge  of 
a primary  beam  particle.  Since  different  materials  have  different  I and 
de/dx  , the  volumetric  ionization  rate  for  a mixture  of  several  species  is 

-des/dx \ 

/ (1-2) 

where  S'  is  the  fractional  concentration  of  species  S 

Secondary  electrons  produced  in  this  manner  emerge  from  the  reaction 
with  a variety  of  energies  thus  creating  a secondary  source  spectrum  of 
electrons  which  depend  on  the  energy,  c , of  the  primary  electrons.  If 
one  defines  a differential  ionization  cross  section  for  a primary  particle 
with  energy  e to  produce  a secondary  electron  with  energy  T as  S(e,T) 
(units  = cm?/eV)  then  the  total  ionization  cross  section  is 

Oi(e)  = S(e,T)dT  (1_3) 

where  Tm  is  the  maximum  kinetic  energy  of  an  ejected  electron.  Tf  one 
assumes  that  the  fastest  outgoing  electron  from  the  reaction  is  the  primary, 
then  Tm  = (e  - I)/2  where  I is  the  ionization  potential  of  the  neutral 
species.  Therefore,  SQ  can  be  rewritten  as 

s°  3 ^ ns  °is 

q s (1-4) 

where  ng  is  the  number  density  of  species  s , and  the  secondary  electron 


^o  3eb  ^ ^s  / 

— ( 


source  spectrum  can  be  expressed  as 


Number  of  electrons  \ 
per  unit  volume  per 
unit  time  produced 
with  energy  between  > =■  i>(T)dt  = ^eb  E nsSs(e,T)dT 

T and  T + dt  by  j q 

primary  electrons 
with  energy  e 


J 


(1-5) 

From  equation  1-2  and  1-4  one  should  also  note  the  relationship  between 
the  ionization  cross  section  and  the  energy  lost  per  primary  electron  per 
unit  length  of  penetration 

Ofs  (e)  = 1 /-drs/dx\ 

”s\  ~*B~)  (1-6, 

This  formulation  has  been  structured  to  determine  the  volumetric  production 
rate  of  secondary  electrons  due  to  a mono-energetic  electron  beam. 

A second  way  of  approaching  the  problem  is  to  analyze  the  complete 
energy  degradation  of  a primary  electron  with  energy  e . Let  the  total 
number  of  secondary  electrons  with  energy  T produced  by  a primary  electron 
of  energy  e be  defined  as 

N (T,  = /*  , f(e',T)de' 

e emrn  (1-7) 

where  f(e’,T)  is  a function  that  involves  the  probability  of  producing  a 
secondary  electron  with  energy  T as  well  as  the  probability  that  the 
primary  electron  has  energy  c'  . The  minimum  energy  the  primary  electron 
may  have  and  still  produce  a secondary  electron  of  energy  T is 

e . = 2T  + I (1-8) 

min 


if  the  fastest  outgoing  electron  in  the  reaction  is  assumed  to  be  the 
primary.  f(E#fT)  may  be  expressed  as 

fte'.Tjde*  = N ) (1-9) 

where  N is  the  number  of  scatters  encountered,  P(ef,T)  is  the  probability 
of  producing  a secondary  electron  with  energy  T,  and  wU*)  is  the 
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probability  that  the  primary  electron  has  energy  c' 
a)(e')  may  be  expressed  as 


P (e' ,T)  and 


P(e'  ,T)  =•  S (e'  ,T) 
unit  area 


(1-10) 


w = dx 


where  R is  the  range  of  the  primary  electron  expressed  as 


(1-11) 


R - fc 


1 de' 
de ' /dx 


(1-12) 


and  dx  is  the  fraction  of  the  range  where  the  primary  electron  has  energy 
e'  . Thus,  equation  1-9  becomes 


f(e',T)  = ft  of  scatters  j S(e',T)  de'  1_ 
per  unit  area  / (dtc'/dx)  R 


(1-13) 


where  the  number  of  scatters  per  unit  area  equals  Rns  where  ns  is  the 
density  of  scatters. 

Equation  1-7  then  becomes 


Ne  = / . S(e',T)  1_ de‘ 

min  1-1  de' \ 


-1  dej_\ 
ns  dx  j 


(1-14) 


The  # of  secondary 
electrons  with  energy 
T produced  by  a pri- 
mary electron  of  energy 
e'  in  length  dx 


= dNe  = ns  S(e' , T) 


(1-15) 


If  the  incident  flux  of  primary  electrons  is  F = jeb/1^  then  equation 


1-5  may  be  written  as 

£> (T)  = jeb  ns  S(e,T) 


(1-16) 


Secondary  electrons  produced  in  this  manner  are  then  free  to  collide 


with  neutral  species  causing  electronic  and  vibrational  excitation  and 
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additional  ionization.  To  examine  the  excitation  produced  by  all  genera- 
tions of  electrons,  the  loss  function  may  be  introduced 


He) 


-1  d£ 
ng  dx 


(1-17) 

which  expresses  the  energy  lost  per  unit  length  of  penetration  per  particle 
per  unit  volume.  The  loss  function  may  alternately  be  written  as 

(e-I) 


L(e)  = E u>.a  . (e)  + E 
i J J i 


We>  + 


/. 


2 TSi(e,T)dT 


(1-18) 


where 


. th 


uk  = excitation  energy  of  j excited  state 
0^  = cross  section  for  excitation  of  state  j 
1^  = ionization  potential  for  state  i 

= cross  section  for  ionization  of  state  i 
Therefore,  the  total  number  of  excitations  to  state  j by  all  generations 
of  electrons  (assuming  all  excitations  take  place  from  the  ground  state) 
becomes 


(e-  I) 

Ja(e)  = /e.°J(e,)  df- * + r 2 


"j 


L(e') 


/. 


Jj (T)n(e,T)dT 


(1-19) 

where  the  first  term  on  the  right  of  equation  1-19  expresses  the  number  of 
excitations  due  to  the  primary  electron  alone  and  the  second  term  expresses 
the  number  of  excitations  resulting  from  secondary  electrons. 

The  reaction  efficiency  can  then  be  defined  as  the  fraction  of  the 
incident  energy  e going  into  exciting  state  j as 


Vi 


' e (1-20) 

The  number  of  electron  volts  required  to  produce  an  electron  pair  is 


1 * ui  = 


n.  Jj(e) 


(1-21) 
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where  I is  the  effective  ionization  potential  discussed  earlier.  Since 
all  inelastic  collisions  do  not  cause  ionization,  np  < 1 and,  in  gen- 
eral, 1 > I, 

The  various  source  terms  can  now  be  related 


So  3 3eb  / -de/dx 


So  3 3eb  ns0i(e) 


(Equation  1-1) 


(Equation  1-4) 


jeb  ns 

q i 


(1-22) 


The  total  rate  of  electron  production  by  an  electron  beam  can  be  found 
by  integrating  SQ  over  the  electron  range  or 


# of  electrons  produced 
per  cm2  of  beam  per  unit 
time  by  beam  electrons 


■*>T+I  (-de'/dx) 


(1-23) 


Sample  loss  functions  and  efficiencies  for  argon  and  mercury  may  be 
found  in  Appendix  A. 

One  should  note  that  in  performing  Boltzmann  calculations  to  determine 
the  excitation  rate  of  a given  species  due  to  an  electron  beam,  either 
equation  1-16  or  1-19  could  be  used.  If,  however,  a dc  field  is  applied 
to  the  plasma,  the  secondary  electrons  are  accelerated  after  being  produced 
and  thus  may  be  capable  of  more  excitation  than  predicted  by  equation  1-19. 
Thus,  in  order  to  compute  the  amount  of  excitation  due  to  energy  supplied 
by  both  the  electron  beam  and  the  applied  field,  it  is  more  appropriate  to 
make  use  of  the  secondary  source  spectrum  of  equation  1-16  to  monitor  the 


time  evolution  of  the  electron  distribution  function. 
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Pi f forenti al  Ionization  Cross  Section 

In  the  last  section,  the  importance  of  the  secondary  source  spect rum 
in  electron  beam  related  plasma  analysis  was  stressed.  However,  while  the 
fate  of  the  primary  or  incident  electrons  in  the  collision  process  is 
fairly  well  understood,  a great  deal  of  uncertainty  surrounds  the  dynamics 
of  secondary  electrons  being  produced. 

Through  the  use  of  scattering  theory  and  some  standard  approximations 
for  many-particle  wave functions,  matrix  elements,  generalized  oscillator 
strengths  and  Bet  he  surfaces,  various  people  have  attempted  to  develop 
analytic  approximations  for  differential  ionization  cross  sections.  "The 
difficulty  here,  however,  is  that  the  results  often  depend  sensitively  on 
the  scattering  approximations  made,  the  assumed  coupling  schemes,  and  the 
atomic  wavefunctions  and  it  is  not  at  all  obvious  what  the  reliability  is 
in  general."  (Ref  6) 

For  this  reason,  and  due  to  the  availability  of  an  ever-increasing 
amount  of  experimental  data,  a number  of  people  have  turned  to  a phenomeno- 
logical approach  in  which  an  analytic  representation  of  t lie  secondary 
source  spectrum  is  developed  by  using  a little  insight  and  a lot  of  trial 
and  error  to  develop  a function  which  agrees  with  available  experimental 
data. 

One  functional  form,  developed  by  Green  and  Sawada  (Ref  3) , shows 
a great  deal  of  promise  when  compared  to  experimental  data  obtained  by 
Opal  et  a_l  (Ref  7)  and  later  used  by  Peterson  and  Allen  (Ref  6)  in  analysis 
of  argon,  and  Bass,  Berg,  and  Green  (Ref  4)  in  analysis  of  mercury. 

This  differential  ionization  cross  section  is  defined  as 

r r>  f 

S(e,T)  =*  A(e)[(T  - ToP  + r7]  (1-24) 
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where  e is  the  energy  of  the  incident  electron;  T is  the  energy  of 
the  secondary  electron  being  produced;  and  A , T0  ' and  ^ are  func_ 
tions  of  e . 

One  strong  point  in  using  this  Lorentzian  form  is  that  when  < = 1, 
it  can  be  easily  integrated  over  all  energies  of  the  secondary  electrons 
to  obtain  a total  ionization  cross  section 


°i(E)  = ;om  s<£'T>dT  = ArLtan"1(Tm  ~ T0>/r  + tan_1(To/r)J  (1-25) 

where,  as  in  the  previous  section, 

Tm  = 1/2  (e  “ (1-26) 

is  the  maximum  energy  a secondary  electron  may  have  if  the  fastest  outgoing 
electron  is  assumed  to  be  the  primary.  Due  to  the  confusion  between  pri- 
mary and  secondary  electrons  at  T > T , attempting  to  examine  the  be- 

m 

havior  of  s(e,T)  for  T > T should  be  viewed  with  considerable  caution. 

m 

A reasonable  form  for  the  resonance  parameter,  Tq  , was  found  by 
Green  and  Sawada  to  be 


T = T 
o s 


e + T, 


(1-27) 


where  Ta  = 1000,  = 21  , and  Tg  is  adjusted  to  fit  the  gas  being 

considered. 

The  width  parameter,  P , was  chosen  to  be  of  the  form 


r = T £ 

s e + rb 

where  T.  = I , in  most  cases,  and  T 
d s 

considered. 


(1-28) 

is  dependent  on  the  gas  being 
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The  form  for  A(e)  was  found  to  be 


A (e) 


(1-29) 


— IS  2 

where  o = 10  cm  and  K and  J are  dependent  on  the  gas  being  con- 
o 

sidered. 

With  those  forms  for  Tq/  T,  and  A , it  is  easy  to  see  that  for 
large  energies 
a . (c)  * lne 

1 e (1-30) 

which  is  the  traditional  asymptotic  form  for  the  ionization  cross  section 
in  the  Born  Bethe  approximation 

Comparisons  between  differential  ionization  cross  section  points  pro- 
jected by  equation  1-24  and  selected  experimental  data  from  Opal  et  al 
(Ref  7)  may  be  found  in  Appendix  A. 
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— - vwamm 


I I . Boltzmann  Equation 


The  Boltzmann  equation  is  one  of  the  principle  tools  used  in  the 
study  of  collision  dominated  systems.  In  the  case  of  plasma  analysis, 
it  attempts  to  relate  changes  in  the  electron  energy  distribution  func- 
tions to  collisions  between  electrons  and  heavy  particles,  the  effects  of 
applied  fields,  changes  induced  by  external  sources  of  highly  energetic 
electrons,  and  the  effects  of  various  electron  loss  mechanisms. 

Although  the  derivation  of  the  Boltzmann  equation  will  not  be  per- 
formed here,  it  will  be  necessary  to  define  the  various  variables  and 
terms  to  be  used  in  development  of  the  analytic  solution  of  Chapter  IV. 

Following  the  guidelines  of  Rockwood  (Ref  5),  the  flux  divergent 
form  of  the  Boltzmann  equation  can  be  expressed  as 


+ R , (e  + e*  . ) n (e 
si  si 


+ e*  . ) 
si' 


+ 


6(e)/  * 
e*  . 
si 


R . (e)n(e)de 
si 


R .(e)  + R . (e)Ns  + R (e)] 

\ s3  s3  ‘o  si  j 

Ns 


£ (e)-6  (e)N  R (e)n(e) 
s a a 


(2-1) 
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. . w.... 


where 


n (e , t) de 


N 

v/N 


q 


s 


m 


e 

T 


R . 
S3 


(O 


R .(c) 
s] 


R . (c) 
si 


s 


s 


V 


* 

e . 
si 


N 

a 

Ra  (e) 

d 


=»  the  number  of  electrons  per  unit  volume  with  energies 
between  e and  e *■  de(l/cm3) 

= total  gas  number  density  (1/em-*) 

= reaction  rate  constant  for  momentum  transfer  (cmVsec) 

=»  mole  fraction  of  species  s 
=»  electron  mass 
= electron  charge 
= gas  temperature 

= rate  constant  at  which  electrons  with  velocity  u(e) 
produce  excitation  from  the  ground  state  of  species  s 
to  excited  state  j , losing  energy  c*  • (cm3/sec) 

S J 

= rate  constant  at  which  electrons  at  energy  c suffer  super- 
elastic collisions  with  atoms  in  state  N-^  and  gain  energy 

s 

e* • (cm3/sec) 

=*3 

= rate  constant  at  which  electrons  of  energy  c ionize  spe- 
cies s while  losing  energy  e^.  (cm3/sec) 

= density  of  species  s (1/cm3) 

= density  of  species  s in  excited  state  j (1/cm3) 

= average  collision  frequency  between  electrons  and  heavy 

particles  for  transferring  energy  (v  = — u) (1/sec) 

M 

= energy  required  to  excite  species  s from  the  ground 
state  to  excited  state  j 

= energy  required  to  ionize  species  s from  the  ground  state 
= number  density  of  dissociative  attachment  species 
= rate  at  which  electrons  with  energy  e react  with  species  a, 
dissociate  a,  and  attach  themselves  to  one  of  the  by-products 
of  the  reaction 
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excited-excited  state  transitions  is  relatively  straightforward  but  the 
notation  is  cumbersome;  hence,  they  were  omitted  for  reasons  of  clarity. 
Equation  2-1  also  assumes  the  number  density  of  species  s in  the  ground 
state  approximately  equals  the  total  gas  density  for  species  s . 

A further  explanation  of  the  various  terms  may  be  helpful. 


Term  1 describes  charges  in  the  distribution  function  due  to  the  accelera- 
tion of  electrons  by  an  applied  dc  field.  n/2e  controls  the  electron  flux 
up  the  energy  axis  and  3n/3e  controls  the  diffusion  of  electrons  as 
they  move  up  the  axis. 


Term  2: 

L 

V 

n ( S ' A 

| - KTe  3 n 

3e  1 

i 

\2  J 

3cJ 

Term  2 describes  changes  due  to  momentum  transfer  occuring  when  electrons 

experience  elastic  collisions  with  heavy  particles  whose  temperature  is  T. 

Term  3:  R_ . (c  + e*  )n(e  + c*  )N° 

SJ  sj  sj  s 

* 

Term  3 describes  gains  at  energy  e due  to  electrons  with  energy  £ + e 

sj 

colliding  with  species  s in  the  ground  state,  exciting  s to  state  j 

while  losing  energy  e* , . 

sg 

if  A -i 

Term  4:  R . (e  - e )n(e  - e ) NJ 
S3  S3  sj  s 

* 

Term  4 describes  gains  at  energy  e due  to  electrons  with  energy  £ - £ . 

SD 

undergoing  superelastic  collisions  with  species  s in  excited  state  j 
and  gaining  energy  e*  while  de-exciting  s to  the  ground  state. 
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Term  5;  R . (e  + e*.)n(e  + c*  )N° 

si  si  si  s 

* 

Term  5 describes  gains  at  energy  e due  to  electrons  with  energy  e + e 

si 

colliding  with  species  s in  the  ground  state,  ionizing  s and  losing 

* 

energy  e . . 

si 

Term  6:  6(e)N°  r .(e) n(e)de 

s e . si 
si 

Term  6 represents  gains  due  to  electrons  with  energies  greater  than  the 

ionization  threshold  energy  for  species  s colliding  with  species  s in 

the  ground  state,  ionizing  s , and  creating  secondary  electrons.  The 

6(e)  in  front  of  this  term  signifies  that  all  secondary  electrons  created 

in  this  manner  are  assumed  to  emerge  with  zero  energy. 

Term  7:  R„.(e)n(e)N° 

* J S 

Term  7 describes  losses  at  energy  e due  to  eiecLLUns  wx-th- 

colliding  with  species  s in  the  ground  state,  exciting  s 

* 

while  losing  energy  egj  . 

Term  8:  R . (e)n(e)N^ 

S3  s 

Term  8 describes  losses  at  energy  e due  to  electrons  with 

colliding  with  species  s in  state  j and  gaining  energy 

de-exciting  s to  the  ground  state. 

Term  9:  R . (e)n(e)N° 

si  s 

Term  9 describes  losses  at  energy  e due  to  electrons  with 
colliding  with  species  s in  the  ground  state,  ionizing  s 
energy  egi  . 

Term  10;  £g(e) 

Term  10  represents  gains  due  to  an  external  flux  of  electrons  colliding 
with  species  s in  the  ground  state,  ionizing  s , and  producing  second- 
ary electrons  with  energy  e . If  the  energy  of  the  incident  electrons 


to  state  j 

energy  e 

e*  . while 
S3 

energy  c 
, and  losing 
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is  sufficiently  large  (.3MeV),  the  change  in  energy  of  the  incident  elec- 
trons can  be  neglected . 

Terra  11:  5(l)N  R (t)n(c) 

— a a 

Term  11  represents  losses  due  to  electrons  with  energy  e reacting  with 
species  a , dissociating  a , and  attaching  themselves  to  one  of  the  by- 
products of  the  reaction.  The  6(e)  symbol  signifies  that  only  electrons 
with  zero  energy  may  be  attached  in  this  manner. 

If  one  limits  his  investigation  to  a quasi  steady  state  system,  time 
dependence  can  be  eliminated.  If  one  assumes  the  density  of  excited  atoms 
is  much  less  than  the  density  of  ground  state  atoms,  superelastic  col- 
lisions may  also  be  ignored.  And,  finally,  if  one  can  assume  that  the 
effects  of  elastic  collisions  are  much  smaller  than  the  effects  of  inelas- 
“tic  collisions, — *■'  ..oiiiginiis  may  also  be  ignored.  With  these 


restrictions,  Rockwood's  equation  ( 2— 1 ) reduces  to 


-2U*  (K/N) 2 3 


3m 


+ E . N 
S,  3 s 


9c 


(" 

lv/N  \2 


2c  dc 


R . (e  + e ,)n(e  + e*  ) + R (c  + t*  )n(e  + e*  ) 
sj  sj  sj  si  si  si 


- 6(e)/  R ,(e)n(e)de  - R ,(c)n(e)-R  .(e)n(e) 
e*  si  sj  sx 

si 


+ E S (e)  - 6(e)N  R (e)n(e)  =*  0 
s s a a 


(2-2) 


further  simplification  can  be  accomplished  by  redefining  the  distribution 
function  from 


to 


/ n(e)de  =*  n » total  electron  number  density 
o © 


/“  e^ftcjde  =■  1 
o 


(2-3) 


(2-4) 
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(2-5) 


where  n(e)  = nee*5f(e)  and  defining 


Q(e)  = Eq  o (e) 

S S •> 


N°  = Nq 


R .(e)  = — | o . (e)  e*5 

s]  (ml  S3 


Rsi<C)  = 'i 


(2-6) 

(2-7) 

(2-8) 


Equation  2-2  ther  becomes 

e2 (K/N) 2 3_  T e 3f(e)l 

3 9e  l_Q(e)  3e  J 

+ s*j  qs  [°sj(£  + EIj)(e  + Elj)f(e  + Esj) 

♦ o . (e  + c .)  (c  + e ,)f(e  + e .) 
si  SI  si  SI 


+ 6(c)/  * a . (e)ef(e)de 
Lsi  si 

-a  . (e)  t f (e)  - o (c) cf  (c ) 
S3  si 


+ IA!iL_  . 

(2/m) ^Nn 


6(e)fa°a(u)uf (u)  = 0 


(2-9) 


where 


= mole  fraction  of  species  a 

= excitation  cross  section  for  species  s to  state  j 
= ionization  cross  section  for  species  s 
o = attachment  cross  section  for  species  a 

d. 

And,  finally,  converting  energy  to  electron  volts  (u  = e/e)  , equation 
2-9  becomes 


'sj 


si 
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(E/N) 2 9 u 9f (u) 

3 2u  Q(u)  9(u) 


+ 


* * * 

E q o . (u  + u .)  (u  + u ,)f  (u  + u .) 
s,j  s sd  sj  S3  S3 


+ 


* * * 

o . (u  + u , ) (u  + u , ) f (u  + u . ) 

SI  si  SI  SI 


+■  <5(u )/  t o . (u)  (u)  f (u)du 

U . SI 
SI 

- a . (u)  (u)  f (u)  - o . (u)  (u)f  (u) 
s3  si 

+ cE£  (u)  - 6(u)Y  o (u)(u)f(u)  = 0 (2-10) 

s s a a 

where 


/2V*Nn  (5.93xl07)Nne 


(2-11) 


This  equation  forms  the  basis  for  the  analytic  work  performed  in  Chapter  IV. 


95 


1 


r 


III.  Numerical  Solution 

In  analysis  of  all  but  the  simplest  collision  dominated  phenomena, 
the  Boltzmann  equation  adequately  expresses  the  various  processes  involved, 
but  proves  much  too  complicated  to  be  solved  analytically.  For  this  rea- 
son, most  groups  involved  in  plasma  related  research  have  turned  to 
Boltzmann  computer  codes  to  analyze  the  time  evolution  of  the  electron 
distribution  function  and  how  it  is  affected  by  variations  in  species 
concentrations,  cross  section  approximations,  applied  fields,  electron 
beams,  etc.  One  such  code  was  developed  in  conjunction  with  this  and 
three  related  theses  by  William  Ercoline,  Thomas  Gist,  and  Harold 
Hastings.  Various  aspects  of  the  Numerical  Gas  Breakdown  (NGB)  program 


will  be  discussed  in  this  chapter.  Details  and  program  documentation  may 
be  found  in  a Technical  Report  (Ref  11)  pending  publication. 

P rogram  Logic 

NGB  considers  the  time  rate  of  change  of  the  electron  distribution 
function  due  to  seven  processes 


dt 


= dn (e ) 

\ + dn(e) 

| 

dn  (e  ) 

dt 

/ dt 

1 

dt 

Non  E-Beam 

Direct 

Electron- 

and  Indirect 

Electron 

Electron 

E-Beam 

Ionization 

Beam 

Ionization 

Collisions 

+ dn(e) 

\ + dn  (e ) ^ 

x dn (e ) 

\ + dn(e)\ 

dt 

/ dt  J 

dt 

I dt  | 

Applied 
dc  Field 


Inelastic 

Excitation 


Dissoci- 

ative 

Attachment 


Elastic 

Collisions 

Between 

Electrons 

and  Heavy 

Particles 

(3-1) 
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The  first  teim  represents  an  electron  source  due  to  free  elections 
with  energy  greater  than  the  ionization  threshold  energy  colliding  with 
heavy  particles  and  producing  secondary  electrons.  All  secondary  elec- 
trons produced  in  this  manner  are  assumed  to  emerge  with  zero  energy. 

The  second  term  represents  an  electron  source  due  to  electrons  from 
mono-energetic  electron  beam  colliding  with  heavy  particles  and  producing 
secondary  electrons.  Secondary  electrons  produced  in  this  manner  emerge 
with  a variety  of  energies. 

Term  three  represents  dispersions  in  the  electron  distribution  func- 
tion due  to  collisions  between  free  electrons. 

Term  four  describes  changes  in  the  distribution  function  resulting 
from  the  acceleration  of  free  electrons  by  an  applied  dc  field. 

Term  five  describes  changes  in  the  distribution  function  due  to  in- 
elastic collisions  between  electrons  and  heavy  particles. 

Ti'ira  six  describes  electron  losses  due  to  electrons  of  zero  energy 
reacting  with  a small  concentration  of  a quenching  species  resulting  in 
the  dissociation  of  the  heavy  particle  and  t lie  attachment  of  the  free 
electron. 

The  final  term  represents  changes  in  the  electron  distribution  func- 
tion due  to  elastic  collisions  between  free  electrons  and  heavy  particles. 

Rach  term  is  computed  independently  of  the  others.  In  this  manner, 
it  is  possible  for  any  combination  of  the  seven  processes  to  be  considered, 
and  the  relative  effect  of  each  term  on  the  distribution  function  can  be 
monitored. 

Method 

As  in  Chapter  II,  the  basic  approach  for  solving  the  dc  Boltzmann 
equation  follows  the  formulation  of  Rockwood  (Ref  5).  The  time  rate  of 
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change  of  the  number  density  n(e,t)de  of  electrons  with  energy  between 
e and  £ + de  in  a mixture  of  gases  of  total  number  density  N is  deter- 
mined by 


3n  - - 3Jf  - 3J„,  + Z.  N°  / R_,(e  + £*.)n(e  + e*.) 
i el  s,]  sis]  sj  sj 

3t  3e  3e  ' 


♦ R •(£  - £*.)n(e  - e*.)Ns  + R (e  + e*  )n(e  + £*  .) 

s J S]  S]  —57  si  si  si 


+ <5  (e)  / * R . (e)  n (e)  de  - [Rci(e)  + R . (e)  + R„i  (e)  ]n  (e)|  (3-2) 

e . si  SJ  sj  si  i 

si  1 

The  first  term  on  the  right  hand  side  of  equation  3-2  is  the  flux  of 

electrons  in  energy  space  driven  by  an  applied  dc  field 


, 2N2(E/N)2 

Jf  = —2 , 

3m(v/N) 


In  - 3n  \ 

l 2T  le  J 


(3-3) 


■ ■ i?r 


| qsos(e) 


where 


qs  = the  mole  fraction  for  species  s 

ag(e)  = the  momentum  transfer  cross  section  for  species  s 

m = electron  mass 

e - electron  charge 

The  second  term  on  the  right  hand  side  of  equation  3-2  is  the  flux 
of  electrons  along  the  energy  axis  driven  by  elastic  collisions 

Jel  • " “(  r ■ e)  - KTe  If]  (3- 


I—)1*  Z qs°s 

/ S Me 


(3-6) 
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gas  temperature 


The  final  term  represents  the  various  inelastic  and  superelastic 


collisions  which  yive  rise  to  non-local  interactions  in  energy  space.  The 


quantity  R -(e)  =»  o,.(e)u(e)  is  the  rate  at  which  electrons  with  velocity 


u(e)  produce  excitation  from  species  s to  excited  state  j while  losing 


Reverse  rates  are  computed  using  the  principle  of 


energy  c 


is  the  rate  at  which  electrons 


ie)  ionize  species  s while  losing  energy  e 


term  multiplied  by  6(e)  indicates  that  all  secondary  electrons  produced 


in  this  manner  emerge  with  zero  energy 


Various  modifications  were  made  to  Rockwood's  basic  formulation  in 


obtaining  the  NGB  Boltzmann  code.  An  electron  source  term  J?(e)  due  to 


an  electron  beam  and  a sink  term  due  to  dissociative  attachment  were  added 


The  source  term  was  computed  following  the  guidelines  discussed  in 


Chapter  I where  the  secondary  electron  source  spectrum  produced  by  a mono 


energetic  electron  beam  is  defined  as 


electron  beam  current  density 
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N = total  gas  density 

Ss(e,T)  = differential  ionization  cross  section  for  species  s 

£p  = energy  of  incident  electrons 

The  dissociative  attacliment  term  was  computed  by 


dn(c ) 
dt 


Attachment 


= 6 (e)NH1  n(e)R,  (e) 

a a 


where  R (e)  = o (e)u(e)  is  the  rate  at  which  electrons  react  with  a 

a a 

molecular  species  resulting  in  the  dissociation  of  the  molecule  and  the 
attachment  of  a free  electron,  is  the  mole  fraction  of  the  quenching 

species.  The  6(e)  indicates  that  electron  attacliment  only  occurs  for 
electrons  with  zero  energy. 

Although  equation  3-2  represents  only  ground  state  to  excited  state 
and  excited  state  to  ground  state  transitions,  the  equation  was  modified 
to  include  excited  state  to  excited  state  transitions. 

Equation  3-2,  with  the  modifications  just  mentioned,  was  converted 
into  a set  of  K coupled  ordinary  differential  equations  by  finite  dif- 
ferencing the  electron  energy  axis  into  K cells  of  width  Ac  and  rear- 
ranging them  such  that  .cLrLi£)  due  to  each  of  the  seven  processes  described 

dt 

in  equation  3-1  were  isolated.  The  manner  in  which  this  was  done  is  des- 
cribed in  Reference  11. 

Contributions  due  to  each  of  the  various  terms  could  be  summed  to 

obtain  a total  which  could  be  passed  to  the  integrator  to  compute 

dt 

the  time  evolution  of  the  electron  distribution  function. 

The  integrator,  DVOGER,  used  by  NGB  was  obtained  from  the  International 
Mathematical  and  Statistical  Libraries,  Inc.  (Ref  8)  and  is  based  on  a 


sixth  order  GEAR  method  (Ref  9) . Along  with  the  distribution  function, 
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NGB  calculates  a number  of  related  parameters  useful  in  analyzing  plasma 
phenomena . 


Pumpi n<j_  Rate  Constants 

Excitation  pumping  rate  constants  are  the  average  rates  at  which 

electrons  produce  excitation  from  species  s to  excited  state  j while 

losing  energy  e*  . . These  pumping  rates  are  computed  in  the  following 
s3 


manner 


<R  .>  =*  — /"  a • (e)u(e)n(e)de 
SJ  n w 


(3-10) 


where  ng  is  the  average  electron  density.  Reverse  pumping  rate  constants 
are  computed  using  the  principle  of  microreversibility 


<R  .> 

S3 


— r 

ne  ° 


e + e . 
sj 


aSj(e  + cs^)u(e)n(c)de 


(3-11) 


Similarly,  ionization  pumping  rate  constants  are  the  average  rate  at 
which  electrons  ionize  species  s while  losing  energy  eV  and  are  com- 


puted by 


<R  • > = _i.  /”  a_<  (c)u(e)n(e)dc 


si 


ne  ° 


si 


(3-12) 


Energy  Absorption  and  Depletion  Rates 

Energy  absorption  and  depletion  rates  are  the  rates  at  which  energy 
is  either  gained  or  lost  from  the  free  electron  gas  due  to  each  of  the 
seven  processes  described  in  equation  3-1.  NGB  computes  absorption  and 
depletion  rates  for  each  process  by 

EDOT  = / c I ^n_(e)  \ de 
dt 
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(3-13) 


is  the  time  rate  of  change  of  the  electron  distribution 

due  to  each  term  in  equation  3-1.  It  should  be  noted  that  in  true  steady 
state 

7 

kdot  = 0 

T 

T=1 

Drift  Ve locity 

Drift  velocity  is  a parameter  used  to  compare  numerical  predictions 
with  experimental  data  and  is  defined  as 

u = -1—  (KDOT) 
a neE  Due  to  Applied 

dc  Field  (3-14) 

where  KDOT  is  the  rate  at  which  electrons  gain  energy  due  to  an  applied 
dc  field,  E is  the  electric  field  strength,  and  n^  is  the  electron 
density. 


where 


m. 


Reaction  Efficiencies 

Reaction  efficiencies  are  the  fraction  of  the  total  input  energy  which 
supports  each  of  the  inelastic  channels.  Excitation  efficiencies  are 
computed  by  NGB  in  the  following  manner 


v 


sj 


e n N <R  > 
sj  e s sj 


(EDOT) 

Due  to  Applied 
dc  Field 


+ (EDOT) 


Due  to 
E-Beam 


-1 


where 


(3-15) 


esj  = excitation  threshold  energy  for  exciting  species 

s to  state  j 

n^  = number  density  of  electrons 

= number  density  of  species  s 
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<Rgj>  = average  rate  for  exciting  species  s to  state  j 

(EDOT)  + (EDOT)  = the  rate  at  which  energy  is  supplied  to  the  plasma 
by  the  applied  dc  field  and  the  electron  beam 
Similarly/  efficiencies  for  ionization  processes  are  defined  as 


v • = e . n N <R  . > 
si  si  e s si 


(EDOT)  + (EDOT) 

Due  to  Applied  Due  to 

dc  Field  E-Beam 


-1 


(3-16) 


where 

= ionization  threshold  energy  for  species  s 
= average  rate  for  ionizing  species  s 
Again,  one  should  note  that  in  steady  state,  energy  balance  requires 


si 


<R  .> 
si 


I (v  + v .)  = 1 
sj  s]  si 


(3-17) 


Rel i abi lity  of  NG B Numerical  Solution 

In  an  attempt  to  verify  the  accuracy  of  results  being  generated  by 
NGB,  a cursory  comparison  was  made  between  drift  velocities  generated  by 
Rockwood  (Ref  5)  and  those  produced  by  NGB  using  similar  input  parameters. 


Figure  3-1  illustrates  the  results. 


E/N  fv  cm*) 

Figure  3-1 

Comparison  of  computed  values  (solid  curve) 
and  experimental  data  from  Refs.  6 and  13-15  for  drift 
velocity  v t (cm/sec)  and  characteristic  energy  (eVj 
as  a function  of  E/N  (Vcm*). , 

(Ref  5:2351) 


The  solid  line  designates  the  Rockwood  results  and  the  x's  represent 
NGB  produced  drift  velocities  using  similar  input  parameters.  The  drift 
velocities  obtained  through  NGB  appear  to  be  consistent  with  those  obtained 
by  Rockwood.  This  correlation  tends  to  lend  credence  to  other  NGB  results. 
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I: 


IV.  Analytic  Solution 


Although  computer  codes  are  capable  of  numerically  solving  the 
Boltzmann  equation,  only  a single  solution  is  produced  for  each  set  of 
physical  parameters  specified.  Hence,  it  is  often  difficult  to  determine 
the  effect  of  each  parameter  on  the  distribution  function. 

Due  to  the  increasing  popularity  of  electron  beam  sustained  lasers, 
beam  effects  on  the  electron  distribution  function  have  gained  considerable 
attention.  To  obtain  some  insight  into  the  problem,  a simple  analytic 
model  was  derived  in  which  emphasis  was  placed  on  the  effects  of  electron 
beam  currents  and  applied  dc  fields  on  the  electron  distribution  function. 

r> f ^rnuhinn  and  analysis  of  a number  of  the  approximations 

required  are  contained  in  Appendices  C and  D. 

Assumptions 

A number  of  simplifying  assumptions  were  made  in  the  analytic  solu- 
tion. Among  these  were 

1)  Steady  state 

2)  Single  species  which  exists  only  in  the  ground  and  singly 
ionized  state 

3)  Secondary  electrons  resulting  from  non-electron-beam  collisions 
emerge  with  zero  energy 

4)  Momentum  transfer  cross  sections  are  independent  of  the  electron 
energy 

5)  Electron  attachment  occurs  only  from  electrons  with  zero  energy 

With  these  assumptions,  the  Boltzmann  equation  discussed  in  Chapter 


II  (equation  2-11)  becomes 


(E^V-  f-  U + c£(u)  + o.(u  + I)  (U  + I)f  (u  + I) 

3Q  du  (_  du  J 1 

- a.  (u)(u)f(u)  - 6 (u) ¥ a (u) (u) f (u)  + 6(u)/r  a . (u) (u) f (u) du  = 0 

1 a a XX 

(4-1) 

where  I is  the  ionization  threshold  energy. 

Since  the  significance  of  many  of  the  above  terms  varies  in  different 
regions  of  the  energy  axis,  solutions  to  the  Boltzmann  equation  above  and 
below  the  ionization  threshold  energy  were  considered  separately  (see 
Figure  4-1) . 


/ 

k 

Region  1:  u > 1 

Region  2:  u < I 

f(u) 

2 

IX. 

u > 

Figure  4-1  Regions  of  Solution 


Constraints 

To  insure  at  least  minimal  validity,  the  following  constraints  were 
imposed  on  the  solution: 

1)  Piecewise  continuity 

f1(I)  = f2(I)  (4-2) 

2)  The  function  and  its  derivatives  must  vanish  at  infinity 

lim  f(u)  = 0 lim  = 0 

u-k»  u-h*>  du  (4-3) 


3)  The  function  must  be  finite  at  u = 0 


(4-4) 


(4-5) 


4)  Electrons  must  be  conserved 

c£(u)  + 6 (u)  /“o  • (u)  (u)  f (u)du  = 6 (u) /“f  o (u) (u) f (u)du 

1 1 Odd 

Reg  ion  1 : (u  >_  I) 

Allowing  for  attachment  only  at  zero  energy  and  assuming  non-electron- 
beam  generated  secondary  electrons  emerge  with  zero  energy,  the  Boltzmann 
equation  for  u greater  than  the  ionization  threshold  becomes 


( E/N ) ^ d_ 
3Q  du 


u df 
du  _ 


+ cS(u)  + o^(u  + I) (u  + I)f (u  + I) 


- oi(u)  (u)f  (u)  = 0 

Approximating  the  inelastic  cross  section  by  (Ref  10) 


(u)  = 0 (1  ~ i) 
° u 


(4-6) 


u o^u)  = oQ(u  - I) 


(4-7) 


(4-8) 


(u  + I)o^(u  + I)  = 0QU 


(4-9) 


where  I is  the  ionization  threshold  energy  and  oq  is  a constant  dependent 

2 

upon  the  species  and  reaction  considered  (approximately  lO^^cm  ) . 

If  one  can  also  assume  that  f (u)  is  not  a rapidly  varying  function 


for  u >_  I then  one  can  further  make  the  approximation 
[ f (u  + I)  - f (u)  ] 3 I 

du 


(4-10) 


Substituting  equations  4-8  through  4-10  into  equation  4-6  yields 


(E/N) ^ d 
3Q  du 


u "1+  or 
_ du  J ° 


d (uf ) . 


-c$  (u) 


(4-11) 


Direct  solution  of  this  second  order  differential  equation  yielded  a form 
which  became  prohibitively  complicated  as  the  solution  in  Region  2 was 
developed.  Therefore,  a modified  solution  was  sought. 
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If 


(E/N) 


3 Q 


d F u d f 1 

du  L du  J 


can  bo  considered  small  compared  to  remaining  terms  in  equation  4-11,  an 
iterative  technique  can  be  employed  to  obtain  an  approximate  solution  in 
which  the  second  order  term  is  ignored.  Using  this  approximate  solution, 
the  second  order  term  can  be  solved  explicitly  and  reinserted  into  equation 
4-11  to  obtain  increasingly  accurate  solutions. 

Using  this  technique,  the  analytic  solution  for  u >_  I becomes 

fx(u)  = cUqSq  1_  + (E/N)2  cUqSq  1 + c" 

,2 


20oI  u" 


3Q  (0oI)2u3 


‘ (4-12) 

where  c''  is  a constant  of  integration  introduced  when  solving  the  dif- 
ferential equation  and  the  secondary  source  spectrum  has  been  approximated 
by 


f S 


S(u)=  / 


o 

2u, 


uS 
o o 


o < u < u^ 
— — o 


V 


2u^ 


(4-13) 

where  SQ  is  the  volumetric  production  rate  of  electrons  due  to  ionization 
from  a mono-energetic  electron  beam 


S = / £(u)du 


(4-14) 


and  u is  the  ionization  threshold  energy  (Reference  Appendix  D) . 
o 

To  satisfy  the  initial  condition  of  normalization 
/“  u^ftujdu  = o 

it  is  clear  that  c' ' must  equal  zero.  Thus,  the  form  of  f (u)  for 
u > I is 
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sm mmmm 


fl(u)  » CUQSQ  1 + (E/N)2  CUQSo  1_ 

3Q  (%D2  u3 


^o1  u2 


(4-15) 


It  should  be  noted  that 


lira  f (u)  = o 

U x» 


and 


lim  = o 


U'X®  du 

as  required  by  the  initial  constraints  on  allowed  forms  for  f (u) . 

Region  2;  (o  u I) 

As  in  Region  1,  if  we  allow  for  attachment  only  at  zero  energy  and 
assume  non-electron-beam  generated  secondary  electrons  emerge  with  zero 
energy  and  recall  that  o^(u)  = o for  u < I , the  Boltzmann  equation 
for  u greater  than  the  ionization  threshold  energy  becomes 


(E/N) 

3Q 


d u df 
du  L du  . 


+ c3(u)  + o^(u  + I) (u  + I) f (u  + I)  = o 


(4-16) 


From  equation  4-9 

(u  + I)o,(u  + I)  = uo 
1 o 

and  from  Region  1 (equation  4-15) 

f (u  + I)  = CUoSo  + (E/N)2  cu0SQ  i 

2ooI  (u+I)2  " 3Q  (0ol)^  (U+I)3 

Substituting  these  values  for  (u  + Do^u  + I)  and  f(u  + I)  into 
equation  4-16  and  solving  for  the  second  order  term  yields 


d j u df 
du  L du  . 


J2£ £(u)  - 3QC  cu0S0 


(E/N) 


- CU 

o o 


(E/N)2  21  (u  + i)7 


OqI 


(u  + I) 


(4-17) 


Again,  unfortunately,  direct  substitution  by  separation  of  variables 
results  in  an  obnoxious  form  for  f (u)  which  appears  to  be  unbounded  at 
u = o . Therefore,  a modified  WKB  was  used  to  express 


(4-18) 


as 

f (u)  = 2/g(u)du  - ug (u)  + const  (4-19) 

(Reference  Appendix  D) . 

Using  this  approach,  equation  4-17  becomes 


f(u)  - - _3Qc — [2/0(u)du  - u8(u) ] 


(E/N)' 


- 12  CUoSo 

2/  u 

du  - u2 

2 - 
(E/N)  21 

(u  + I)2 

(u  + I)2_ 

cu  s 

o o 


21  u 


du 


(u  + I)3  (u  + I): 


+ c 


(4-20) 


Substituting  the  approximation  for  this  secondary  source  spectrum  (equation 
4-13)  into  equation  4-20  and  performing  the  necessary  integrations  yields 
the  solution  to  the  Boltzmann  equation  for  o < u <_  I 
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To  complete  the  solution  for  f (u)  , c' ' ' must  be  determined  and  manipu- 
lations necessary  to  meet  the  initial  noimal ization  requirement  must  be 
performed. 

Apply  Piecewise  Continuity  of  f(u)  to  Compute  c"  ' 

c"'  is  obtained  by  applying  the  piecewise  continuity  boundary 
condition 

fL(I)  = f2  (I)  (4-22) 

Substituting  equation  4-21  for  f^fu)  and  4-15  for  f2(u)  , evaluating 

both  functions  at  I and  solving  for  c" ' yields 


(4-23) 


Normalize 

To  complete  the  solution  for  f (u)  , the  normalization  requirement 
must  be  fulfilled.  An  effort  was  made  to  utilize  the  inherent  normaliza- 
tion in  the  definition  of  f (u) . In  doing  so,  it  was  necessary  to  obtain 

the  steady  state  electron  density,  n , which  appears  in  the  definition 

6 

of  c (equation  2-11)  . 

To  maintain  a steady  state  system,  electron  sources  must  be  balanced 
by  electron  sinks.  Therefore,  the  number  of  electrons  gained  through 
electron  beam  and  field  driven  ionization  must  be  countered  by  electron 
losses  through  attachment  processes.  Thus,  from  equation  4-1 

cyu max  + /°°o  . (u)  (u)  f (u) du  = M1  / o (u)(u)f(u)du  (4-24) 

o oi  a o a 


*/  / / _ cu  S . Al.  2 cu  S 

c = - 3 o o + (E/N)  o o 


8 0QI- 


+ £S a 

(E/N)2  2 


3 Q «oI5 


+ 2uo  ln(2 I)  + 3u0 
I 41 
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where  V is  the  mole  fraction  of  the  species  attaching  electrons,  o 

cl  3 

is  the  attachment  cross  section,  and  u is  the  maximum  energy  of  a 

max 

secondary  electron  produced  via  direct  collision  with  the  electron  beam. 

In  solving  this  equation  for  the  steady  state  electron  density,  n , 
it  was  found  that  the  second  integral  of  equation  4-24  is  infinite  when 
evaluated  between  the  limits  of  zero  and  infinity.  This  is  due,  in  part, 
to  the  approximation  used  for  the  ionization  cross  section  (equation 
4-7) . Where  for  a true  cross  section 

lim  (u)  = o (4-25) 

U K» 

the  approximated  form  becomes 

lim  (j.  (u)  = constant  (4-26) 

u *-°°  1 

Therefore,  a finite  limit,  u*  , was  placed  on  the  upper  limit  of  the 
integral  to  obtain  a finite  value  for  the  field  driven  source  term.  Un- 
fortunately, this  normalization  technique  proved  to  be  very  sensitive  to 

* 

the  value  selected  for  this  artificial  parameter  u . 

Thus,  in  introducing  the  various  approximations  required  to  obtain 
f (u)  its  inherent  normalization  has  been  lost.  As  a consequence  of  this 
loss  of  inherent  normalization,  the  standard  normalization  procedure  was 
employed  in  which  a normalization  constant  defined  by 

/“usf(u)du  = cnom  (4-27) 

was  determined.  The  normalized  f(u)  is  thus  1/c  f (u)  . It  should 

norm 

be  noted  that  this  process,  though  relatively  straightforward,  is  a bit 

tedious.  Evaluating  equation  4-27  yields 

c = /1u1,f,(u)du  + /“u^f.  (u)du  (4-28) 

norm  o * I i 
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Substituting  the  forms  for  f^(u)  and  f (u)  (equations  4-15  and  4-21) 


equation  4-28  becomes 


(4-29) 


where  c" ' is  defined  as  before  (equation  4-23). 


In  obtaining  this  solution  for  c , the  following  approximation 

norm 


was  used 


/Iu1,ln(u  + I)du  5 /V  ln2I  + u - 1 du 
o o — — 

21  2 


(4-30) 


Equations  4-15  and  4-21  along  with  the  normalization  constant  4-29 
form  the  analytic  solution  for  f (u)  for  o < u < 00  . 


"1 
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Analysis  of  Results 


Analytic  Solution 

To  verify  the  accuracy  of  the  analytic  solution  developed  in  Chapter 
IV,  the  Boltzmann  computer  code  (NGB)  was  modified  to  provide  numerical 
solutions  consistent  with  the  various  analytic  assumptions  and  approxi- 
mations made  in  Chapter  IV.  The  Lorentzian  secondary  electron  source 
spectrum  was  replaced  with  the  approximate  form  (equation  4-13) 


£(u)  = 


where  the  volumetric  rate  of  production  of  electrons  due  to  electron  beam- 
neutral  particle  collisions,  SQ  , was  obtained  from  Appendix  D (equation 
D-15)  . 

Argon  was  selected  as  the  test  species  with  an  ionization  threshold 
energy  of  15.7  eV  and  an  energy  independent  momentum  transfer  cross  sec- 
tion of  10“*^cm2.  Ionization  cross  sections  were  replaced  with  the  analytic 
form  (equation  4-7) 

Oi(u)  =o(l-  I) 
u 


where  the  value  for  oQ  was  chosen  to  be  6.9  x 10-1&cm2  (Ref  10).  And 
finally,  excitation  and  inelastic  collision  terms  were  deleted. 

By  varying  the  electron  beam  current  density,  jeb  , and  the  dc 
field  strength,  (E/N)  , a number  of  numerical  solutions  were  obtained 

and  compared  with  their  analytic  counterpart.  Figure  5-1  illustrates 
the  numerical  and  analytic  solutions  obtained  for  a beam  current  density 
of  .001  amp/cm2  and  an  E/N  value  of  10  x 10“17V-cm2  . 


1 
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The  most  immediate  observation  is  that  the  procedure  used  in  Chapter 
IV  has  failed  to  yield  the  desired  normalization.  The  cnorm  factor 
appears  to  be  an  order  of  magnitude  too  large.  Unfortunately,  this  nor- 
malization problem  could  not  be  resolved  in  time  to  be  included  in  this 
thesis.  Hence,  a detailed  analysis  of  the  validity  of  the  analytic  solu- 
tion could  not  be  performed. 


The  approximation  for  the  logrithmic  integral 


f1  u^ln (u  + I)du  = f1  u^T11121  + - i]du 

o o [ 21  2J 

used  in  obtaining  the  normalization  constant  may  partially  be  responsible 
for  the  normalization  problem.  However,  it  does  not  appear  that  this  is 
the  entire  problem  since  the  approximation 

In (u  + I)  = ln2I  + “ A 

21  2 

has  a maximum  error  of  approximately  7%  over  the  range  o <_  u <_  I for 
I = IS. 7 . Thus,  it  appears  that  there  may  be  additional  errors  in  the 

development  of  cnorm  • Due  to  the  considerable  algebraic  manipulations 
required  tc  obtain  cnorm  , a simple  arithmetic  error  cannot  be  completely 
ruled  out. 

Ignoring  the  final  magnitude  of  c and  considering  only  its 

norm 

variable  dependence,  one  can  see  that  the  analytic  solution  is  independent 
of  the  electron  beam  current  density.  Since  f (u)  (equations  4-15,  4-21, 
and  4-23)  depends  linearly  on  SQ  and  the  normalization  factor,  cnorm 
(equation  4-29),  is  also  linear  on  SQ  , f<u^/cnorm  is  indePendent  of 
the  volumetric  production  rate  of  electrons  due  to  the  electron  beam. 

This  independence  is,  in  part,  due  to  the  assumption  used  in  developing 
the  solution  in  the  high  energy  region. 
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In  applying  the  WKB  method,  it  was  assumed  that  field  effects  were 
much  smaller  than  beam  effects.  Since  the  solution  in  Region  1 (u  >_  I) 
was  used  in  developing  the  solution  in  Region  2 (o  <_  u I) , it  is  not 
difficult  to  see  how  this  linear  dependence  on  S0  could  result.  As  a 
consequence  of  the  small  field,  large  beam  assumption,  the  analytic  solu- 
tion is  not  a good  approximation  in  high  energy  regions  for  small  beams. 
Figure  5-2  illustrates  numerical  and  analytic  solutions  for  beam  current 
densities  of  .001,  1,  and  10  amps/cm2  and  a field  strength  of  10xl0-1 7V-cm2 
From  Figure  5-2,  one  can  see  that  as  the  beam  current  density  increases, 
the  analytic  solution  approaches  that  of  the  numerical  solution  at  ener- 
gies above  the  ionization  threshold  energy.  The  analytic  solution,  however 
does  seem  to  be  a reasonable  approximation  to  the  numerical  solution  in 
the  low  region  of  the  energy  axis  where  field  effects  dominate. 

Numerical  Solution 

Using  the  NGB  Boltzmann  code  to  obtain  numerical  solutions,  the 
effects  of  beam  currents  on  a more  complex  system  was  examined.  The  sys- 
tem selected  was  modeled  after  one  used  by  Stephen  Rockwood  in  the  analysis 
of  mercury  (Ref  5) . The  inelastic  processes  considered  were  excitation 
from  the  Hg  ground  state  to  the  63PQ  ^ 2 states,  the  state,  a 

lumped  state  representative  of  higher  electron  levels,  and  ionization 
from  the  ground  state.  Elastic  processes  include  both  electron-electron 
and  electron-heavy  particle  elastic  collisions.  The  dominant  electron 
loss  mechanism  was  assumed  to  come  from  dissociative  attachment  to  Cl2- 
The  initial  electron  distribution  was  assumed  to  be  a Maxwellian  with  an 


average  electron  energy  of  1.75  eV 


(eV) 

Figure  5-2  Numerical  and  Analytic  Solutions  for  Argon 
(E/K  = 10xl0~17V-cm2  and  j . = .001,  1,  10  amp/ cm2) 


Solutions  were  obtained  at  a variety  of  values  of  E/N  and  jgb 
In  all  cases,  however,  problems  developed  in  getting  NGB  to  run  to  a 
steady  state  or  even  a self  reproducing  solution.  After  a closer  exami- 
nation of  the  Boltzmann  equation,  this  failure  to  arrive  at  a self 
reproducing  form  becomes  a bit  more  obvious. 

If  one  assumes  that  the  effects  of  electron-electron  collisions  are 
small  compared  to  other  terms  in  the  Boltzmann  equation  (a  reasonable 
assumption  for  electron  densities  below  1012cm~3)  the  Boltzmann  equation 
can  be  written  in  very  general  terms  as 

. dn  (e) 


dn(e)  _ dn (e) 
dt  dt 

dn  (e) 
+ dt 


+ dn(e> 
dt 


Field 


dt 


Inelastic  Excitation 


Ionization 


Attachment 


dn  (e) 
dt 


Momentum  Transfer 
to  Heavy  Particles 


+ S(e> 


Electron  Beam 


(5-1) 


All  terms  in  this  equation  are  linearly  dependent  on  n (s ) 
source  term  due  to  the  electron  beam. 

To  obtain  a self  reproducing  form 

1 dnx(e)  _ l dn2(c) 
dt  ne,  2 dt 


except  the 


(5-2) 


for  all  e where  n . equals  the  total  number  density  of  electrons  at 

e , i 

time  t,  ne  2 equals  the  total  number  density  of  electrons  at  time 
t + At,  n^(e)de  equals  the  number  density  of  electrons  with  energy  e 
at  time  t , and  n2(e)dc  equals  the  number  of  electrons  with  energy  e 
at  time  t + At.  If  one  defines  ne  2 = Yne  i and  n2^e)  = Ynq(t)  where 
y is  some  constant  depending  on  At  then  equation  5-2  can  be  written  as 


44 


dn^ (e) 

+ 

dn^ (c) 

dt 

Field 

dt 

Inelastic  Excitation 

dn^ (e) 

dn  ^ ( e ) 

dt 

Ionization 

dt 

Momentum  Transfer 

to  Heavy  Particles 


dn  (e) 

t 0(e) 

dt  ~ 

Attachment 

Electron  Ream 

1 

dynx (e) 

dy  n 1 ( e ) 

*ne,l 

dt 

Field 

dt 

Inelastic  Excitation 

■f 

dy  n ^ ( e ) 

dyn^ (e) 

*-  . . 

dt 

Ionization 

dt 

Momentum  Transfer 

to  Heavy  Particles 
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dt 
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Electron  Beam 

Simplifying  yields  the  result 

£?(c)  = yj?(e)  (5-4) 

Equation  5-4  indicates  that  a self  reproducing  form  cannot  be  obtained  unless 
y = 1 , implying  a constant  electron  number  density.  This  seems  to  indi- 

cate that  a steady  state  solution  can  be  obtained  but,  due  to  the  non- 
linearity of  the  electron  beam  term,  a self  reproducing  solution  in  which 
the  electron  density  is  changing,  cannot  be  obtained. 

The  parameter  combinations  used  in  analyzing  beam  effects  unfortunately 
did  not  yield  a steady  state  solution.  Therefore,  in  order  to  compare 

results  for  various  values  of  E/N  and  j , , all  numerical  solutions 

eb 

were  obtained  at  a fixed  point,  10-7  seconds,  in  the  evolution  of  the 
distribution  function. 
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Figure  5-3  illustrates  the  percent  of  input  energy  (supplied  by  the 
dc  field  and  the  electron  beam)  supporting  each  of  the  energy  loss  me- 
chanisms for  a beam  current  density  of  .001  amp/cm^  and  various  values 
of  E/N  . The  results  at  E/N  values  below  3xlO~17V-cm2  should  be 
viewed  somewhat  cautiously.  At  t = 10-7  seconds  into  the  evolution 
of  the  distribution  function  the  average  electron  energy  had  changed  from 

2 

the  initial  value  of  1.75  eV  to  a value  of  . 5 eV  for  j = .001  amp/em  . 

eb 

Various  convergence  indicators  used  by  NGB  to  determine  if  a steady  state 
or  self  reproducing  solution  has  been  reached  indicated  that  considerable 
variations  were  still  occuring  at  t = 10  7 seconds  into  the  evolution 
process  for  small  E/N  values.  These  variations  were  less  pronounced 
for  large  beam  current  densities  at  similar  values  of  E/N  . It  appears 
that  10“7  seconds  was  not  sufficient  time  for  the  initial  distribution 
function  to  adjust  from  an  average  energy  of  1.75  eV  to  an  average 
energy  of  approximately  .5  eV  for  small  field  and  small  beam  current 
densities.  This  accounts  for  much  of  the  unusual  behavior  of  the  effi- 
ciency  results  at  E/N  less  than  3x10  V-cm4  . 

From  Figure  5-3  one  can  see  that  the  dominant  energy  loss  mechanism 
at  low  E/N  values  is  due  to  momentum  transfer  (elastic  collisions) 
between  electrons  and  heavy  particles.  As  the  magnitude  of  the  applied 
field  is  increased,  the  electrons  gain  energy  and  eventually  penetrate 
the  excitation  threshold  energies  for  the  first  and  second  excitation 
levels  (Hg*  and  Hg**)  of  mercury.  At  this  point,  the  percent  of  input 
energy  supporting  these  two  inelastic  channels  begins  to  increase  and 
the  percent  of  input  energy  supporting  elastic  collisions  begins  to  decline. 
However,  as  the  average  electron  energy  increases  as  the  field  strength 
increases,  one  would  expect  the  efficiencies  for  the  third  and  fourth 
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Figure  5-3  Efficiencies  for  Energy  Loss  Processes  as  a 
Function  of  Field  Strength  for  Mercury  (j  . = .001  amp/cm2) 


excited  states  (Hy***  and  Hy****)  and  the  ionization  state  to  also  in- 
crease. But  from  Figure  5-3  one  can  see  that  for  the  Hy***  and  Hy**** 
excited  states  and  the  ionization  state,  the  efficiencies  show  an  initial 
drop  as  E/N  increases  and  then  an  eventual  rise. 

This  drop  is  a direct  consequence  of  the  applied  electron  beam. 

The  beam  maintains  a finite  number  of  electrons  with  energies  greater 
than  the  Hg***,  Hg****,  and  ionization  threshold  energies  independent 
of  the  field  strength  and  tends,  at  these  energies,  to  mask  the  effects 
of  increases  in  the  field  strength.  Thus,  for  small  E/N  values,  as 
the  field  strength  increases  the  amount  of  energy  supporting  higher  exci- 
tation and  ionization  processes  is  primarily  due  to  the  beam  and  remains 
nearly  constant.  As  the  field  grows,  the  total  energy  supplied  to  the 


plasma  grows  while  that  supporting  higher  excitation  and  ionization  remains 
nearly  constant;  thus,  the  percent  of  energy  supporting  these  processes 
decreases.  As  the  field  strength  is  increased  further,  the  field  beyins 
to  accelerate  more  and  more  electrons  to  these  higher  energy  states  until 
eventually  the  field  maintains  a number  of  electrons  in  these  higher  energy 
states  comparable  to  the  number  maintained  by  the  beam.  At  this  i>oint. 


the  efficiency  for  the  higher  excitation  or  ionization  process  beyins 
to  increase  with  further  increases  in  E/N 

Figures  5-4  through  5-11  illustrate  these  results.  Figures  5-4 
through  5-7  show  the  ratio  of  field  effects  to  beam  effects  for  the 
Hg**,  Hg***,  Hg****,  excitation  processes  and  the  ground  state  ionization 
process.  This  ratio  is  defined  as 


U)u(e)dn de 

R dt  Field 

j ,«■>  j,  (c)u  (e)dn(e) 

3 dt  , de 


[ Electron 
Beam 


48 


0.00  20.00  40.00  60.00  80.00  100.00 

E/N 

(lxl0_17v-cm2) 


Figure  5-4  Ratio  of  Field  to  Beam  Effects  for  the  Hg**  Excitation 
Process  as  a Function  of  Field  Strength  (jeb  = .001  amp/cm2) 
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Figure  5-5  Ratio  of  Field  to  Beam  Effects  for 
Process  as  a Function  of  Field  Strength  (jeb 
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Figure  5-6  Ratio  of  Field  to  Beam  Effects  for  the  Hg****  Excitation 
Process  as  a Function  of  Field  Strength  (jeb  = .001  amp/cm2) 


whore  is  the  excitation  or  ionization  threshold  energy  and  is 

the  excitation  or  ionization  cross  section  for  process  j 

These  ratios  are  plotted  as  a function  of  E/N.  Figures  5-8  through 

5-11  illustrate  the  efficiencies  for  the  Hg*‘,  Hg***,  Hg****  excitation 

and  ground  state  ionization  processes  as  a function  of  E/N  for  a 
2 

.001  amp/cm  beam.  As  one  can  see,  the  point  at  which  the  efficiencies 
for  a given  process  begins  to  increase  corresponds  closely  to  the  E/N 
value  at  which  the  field  effects  for  that  process  begin  to  be  comparable 
to  beam  effects  (the  point  at  which  the  field  to  beam  ratio  is  approxi- 
mately equal  to  .1). 

Figure  5-12  shows  the  efficiencies  for  the  Hg****  process  as  a func- 
tion  of  E/N  for  beam  current  densities  of  .001  and  1 amp/cm  . As  the 
beam  current  density  increases,  the  point  at  which  the  field  effects 
become  comparable  to  beam  effects  also  increases.  Thus,  one  can  see 
from  Figure  5-12  that  the  E/N  value  at  which  the  efficiency  for  the 
Hg****  excitation  process  begins  to  increase  is  greater  for  a large  beam 
than  for  a small  beam. 

One  can  also  see  from  Figure  5-12  that  at  large  field  strengths 
(greater  than  70x10  'V-cnr)  the  field  effects  begin  to  dominate  beam 

2 

effects  in  exciting  the  Hg****  level  of  mercury  and  the  .001  and  1 amp/cm 
efficiency  curves  begin  to  merge. 

The  field  to  beam  effect  ratio  also  accounts  for  the  strange  behavior 
of  the  pumping  rate  constants  for  the  Hg**,  Hg***,  Hg****  excitation  pro- 
cesses and  the  ionization  process.  Figure  5-13  through  5-16  illustrate 
the  value  of  various  pumping  rate  constants  as  a function  of  E/N  for 
a .001  amp/cm2  beam.  One  can  see  from  Figures  5-13  through  5-16,  that 
the  rate  constants  for  a given  process  remain  nearly  constant  as  E/N 


53 


Efficiency  (Percent  of  Total  Input  Energy) 


Efficiency  (Percent  of  Total  Input  En- 


Figure  5-9  Efficiency  for  the  Hg***  Excitation  Process  as 
a Function  of  Field  Strength  (j  . = .001  amp/cm2) 
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Figure  5-10  Efficiency  for  the  Hg****  Process  as  a 
Function  of  Field  Strength  (j  b = .001  amp/cm2) 
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Figure  5-12  Efficiency  for  the  Hg****  Excitation  Process  as  a 
Function  of  Field  Strength  (jeb  = .001  and  1 amp/cm2) 
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Figure  5-13  Rate  Constant  for  the  Hg**  Excitation  Process  as  a 
Function  of  Field  Strength  (jeb  = .001  amp/cm2) 
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Figure  5-14  Rate  Constant  for  the  Hg***  Excitation  Process  as  a 
Function  of  Field  Strength  (j  , = .001  amp/cm2) 
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increases  until  the  field  effects  become  comparable  to  beam  effects  above 
the  threshold  energy  for  that  process.  Further  increases  in  E/N  then 
result  in  increases  in  the  rate  constant.  If  the  beam  current  density 
is  sufficiently  large,  the  field  effects  over  a large  range  of  E/N  values 
can  be  neglected  for  processes  with  large  threshold  energies.  One  can 
see  from  Figure  5-17  that  for  a 10  amp/cm2  beam,  E/N  va4lues  ranging 
from  lxlO-^  to  80x10“  -cm  ^ have  little  effect  on  the  magnitude  of  the 

ionization  rate  constant. 

From  the  analysis  of  field  to  beam  effects  in  supporting  various 
excitation  and  ionization  processes,  one  can  conclude  that  analytic  solu- 
tions, such  as  the  one  derived  in  Chapter  IV,  which  assume  that,  at  high 

I 

energies  beam  effects  are  much  greater  than  field  effects  for  large  beams 
(1-10  amp/cm2)  appear  to  be  valid  for  a wide  range  of  field  strengths. 

In  practice  this  implies  that  in  determining  the  optimum  field  strength 
for  exciting  the  upper  level  of  an  electron  beam  sustained  laser,  a wide 
range  of  E/N  values  could  be  used  without  significantly  affecting 
the  production  rate  of  secondary  electrons  if  a reasonably  large  beam 
was  used. 
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VI.  Conclusions  and  Recommendations 

Analytic  Solution 

The  analytic  solution  derived  in  Chapter  IV  displays  some  major  in- 
consistencies when  compared  to  a numerical  solution  developed  for  a similar 
system.  The  analytic  solution  is  too  small  around  zero  energy,  too  large 
at  high  energies  and  does  not  appear  to  be  normalized.  Two  approaches  were 
used  in  an  attempt  to  normalize  the  solution. 

The  first  approach  was  to  use  the  inherent  normalization  in  the  defi- 
nition of  f (u)  . However,  in  obtaining  the  steady  state  electron  density 

required  to  complete  normalization,  it  was  found  that  the  source  of  second- 
ary electrons  due  to  the  applied  dc  field  was  infinite.  This  may  have 
been  a result  of  the  ionization  cross  section  approximation  being  used. 
Therefore,  the  inherent  normalization  in  f (u)  was  lost  in  the  many  ap- 
proximations required  to  obtain  the  analytic  solution. 

The  second  normalization  approach  was  a straightforward  integration 
approach  where  a normalization  constant  was  obtained  by  integrating  u^f (u) 
from  u = 0 to  u = “ . In  doing  so,  some  additional  approximations  were 

required  to  evaluate  the  integrals.  In  introducing  these  additional  approx- 
imations, the  normalization  may  again  have  been  lost;  however,  the  possi- 
bility of  a simple  error  in  integration  has  not  been  completely  ruled  out. 

The  analytic  solution  is  independent  of  the  volumetric  production  rate 
of  electrons  due  to  the  electron  beam.  This  is  a consequence  of  the  "big 
beam  approximation"  used  in  obtaining  the  solution.  In  comparing  the 
analytic  solution  to  a numerical  solution  for  a similar  system,  it  was 
found  that  the  analytic  solution  is  a poor  approximation  of  the  numerical 
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solution  at  high  energies  and  small  beam  current  densities.  However,  as 
the  beam  current  density  increases,  the  analytic  solution  more  closely 
resembles  the  numerical  results. 

Numerical  Solution 

Numerical  solutions  were  obtained  using  a system  of  parameters  modeled 
after  those  used  by  Rockwood  (Ref  5)  in  the  analysis  of  drift  velocities 
and  reaction  efficiencies  for  mercury. 

Reaction  efficiencies  and  pumping  rate  constants  were  obtained  for 
a variety  of  values  of  E/N  and  jel3  . Due  to  the  non-linear  dependence 
of  the  electron  beam  source  term  in  the  Boltzmann  equation,  self  reproducing 
solutions  could  not  be  obtained.  However,  to  compare  results  at  varying 
values  of  E/N  and  jg^  , numerical  solutions  were  obtained  at  a fixed 
point,  10-7  seconds,  in  the  evolution  of  the  electron  distribution  function. 

In  analyzing  the  behavior  of  reaction  efficiencies  for  the  higher 
excited  states  and  first  ionized  state  of  mercury,  the  efficiencies  showed 
an  unexpected  drop  and  then  the  characteristic  rise  as  the  strength  of  the 
applied  dc  field  increased.  This  dip  was  directly  attributed  to  electron 
beam  effects  at  high  energies.  These  beam  effects  at  high  energies  were 
also  responsible  for  the  relative  insensitivity  of  pumping  rate  constants 
for  processes  with  high  threshold  energies  to  variations  in  the  strength 
of  applied  dc  fields.  It  was  found  that  the  pumping  rate  constant  for  a 
given  process  would  remain  fixed  as  E/N  increases  until  the  effects  due 
to  the  applied  field  were  comparable  to  the  effects  of  the  electron  beam 
above  the  threshold  energy  for  that  process.  Only  then  did  increases  in 
E/N  result  in  increases  in  the  pumping  rate  constant.  This  seems  to 
imply  that,  in  determining  the  optimum  field  strength  for  exciting  the 
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upper  level  of  an  electron  beam  sustained  laser,  a wide  range  of  E/N 
values  could  be  used  without  significantly  affecting  the  production  rate 
of  secondary  electrons  if  a reasonably  large  beam  was  used. 


Recommendations 

A number  of  improvements  could  be  made  to  the  analytic  solution. 

To  minimize  normalization  difficulties,  a different  form  of  the  ionization 
cross  section  should  be  used.  Although  the  cross  section  used  in  Chapter 
IV  appears  to  be  a good  approximation  around  the  ionization  threshold 
energy,  it  does  not  drop  off  properly  at  high  energies.  There  are  two 
means  by  which  an  improved  cross  section  could  be  introduced.  Either  a 
completely  different  cross  section  could  be  used  or  a three  region  solu- 
tion could  be  developed  in  which  the  form  of  the  cross  section  in  the 
high  energy  region  could  be  approximated  by  a simple  function  such  as 
a'/u  , where  is  chosen  to  match  cross  section  boundary  conditions. 

The  cross  section  around  threshold  could  be  approximated  as  in  Chapter  IV 

as  a (1  - i ) and  the  cross  section  in  the  low  energy  region  would  be 
u 

zero.  However,  the  solution  developed  in  Chapter  IV  is  closely  tied  to 
the  approximation  used  for  the  ionization  cross  section;  thus,  any  drastic 
changes  to  this  approximation  would  most  likely  produce  a solution  con- 
siderably different  to  the  one  obtained  in  this  thesis. 


1 


67 


w ■ 

S 


Bibliography 


1.  Beyer,  William  H.  Standard  Mathematical  Tables  (Twenty-fourth  Edition) . 

Cleveland:  CRC  Press,  1976. 

2.  Mott,  Nevill  F.  and  H.  S.  W.  Massey.  The  Theory  of  Atomic  Collisions 

(Second  Edition) . Oxford  at  the  Clarendon  Press,  1949. 

3.  Green,  A.  E.  S.  and  T.  Sawada.  "Ionization  Cross  Sections  and 

Secondary  Electron  Distributions,"  Atmospheric  and  Terrestrial 
Physics,  34:1719-1728  (1972). 

4.  Bass,  J.  N.,  R.  A.  Berg,  and  A.  E.  S.  Green.  "Electron  Impact  Cross 

Sections  for  the  Element  Mercury,"  Journal  of  Physics  Bulletin, 

5.  Rockwood,  Stephen  D.  "Elastic  and  Inelastic  Cross  Sections  for 

Electron-Hg  Scattering  from  Hg  Transport  Data,"  Physical  Review  A, 

8 (5) : 2343-2358  (1974). 

6.  Peterson,  J.  R.  and  J.  E.  Allen,  Jr.  "Electron  Impact  Cross  Sections 

for  Argon,"  The  Journal  of  Chemical  Physics,  56  (12) :6068-6076 
(15  June  1972) . 

7.  Opal,  C.  B.,  E.  C.  Beaty  and  W.  K.  Peterson.  JILA  Report  No.  108, 

University  of  Colorado,  Boulder,  26  May  1971. 

8.  IMSL  LIB03-0006V1 . Houston,  Texas:  International  Mathematical 

and  Statistical  Libraries,  Inc.  (July,  1977) . 

9.  Gear,  C.  W.  Numerical  Initial  Value  Problems  in  Ordinary  Differen- 

tial Equations.  Prentice-Hall,  1971. 

10.  Hunter,  A.  M.  II.  Capt,  USAF  (private  communication) . Air  Force 

Institute  of  Technology,  Wright-Patterson  AFB , Ohio,  1978. 

11.  Technical  Report  for  Computer  Program  NGB.  (Publication  pending)  . 

12.  Hunter,  A.  M.  II.  Capt,  USAF.  "Notes  on  Atomic  and  Molecular  Excita- 

tion Through  Collisions  With  Charged  Particles,"  unpublished, 
undated. 

13.  Hunter,  A.  M.  II.  Capt,  USAF.  "Notes  on  Rockwood  Equation,"  unpublished, 

undated. 

14.  Dwight,  Herbert  B.  Tables  of  Integrals  and  Other  Mathematical  Data 

(Revised  Edition).  New  York:  The  Macmillan  Company. 


68 


Appendix  A 


Differential  Ionization  Cross  Sections 

Figures  A-l  through  A-3  illustrate  analytic  fits  obtained  by  Green 
and  Sawada  (Ref  3)  to  selected  experimental  data  compiled  by  Opal  et  al 
(Ref  7) . 

The  analytic  form  for  the  differential  ionization  cross  sections 
used  is  defined  as 


S(c,T)  = A(e) 


(T  - Tq)2  + r2 


where  e is  the  energy  of  the  primary  electron  and  T is  the  energy  of 
the  secondary  electrons  being  produced.  A , T , and  Tq  are  functions 
of  £ and  are  defined  as 


A(e)  = o — in  | e 
° e 7 

i J 


r = r e 
S e + 


T - T 

S a 

e + T. 


— 16  2 

where  oQ  equals  10  cm  and  remaining  parameters  are  adjusted  to  obtain 
the  desired  fit.  Parameters  used  in  obtaining  the  fits  illustrated  in 
Figures  A-l  through  A-3  are  contained  in  Tables  A-l  and  A-2. 


S',  o 


Table  A-l 

A (K).  l'(K)  and  3’0(/i)  lie,  X,  nnd  0,*  f 


E 

A 

r 

T, 

Ho  (/ 

« 24  5 cV) 

r,o 

0 0331 

14  5 

-10  9 

100 

0 0787 

10  1 

-13  8 

200 

0 0270 

13  2 

- 4 30 

300 

0 0175 

1 4-3 

- 2 19 

500 

0 0141 

14-2 

- 3 85 

1000 

0 00530 

10  3 

- 2-57 

2000 

0 00355 

10  2 

- 6-11 

K,  (/ 

==  15  0 eV) 

50 

0 148 

20  0 

— 10  7 

100 

0 177 

12  8 

-0  617 

200 

one 

11-9 

2 54 

300 

0 0SG7 

12  5 

2 40 

600 

0-0G3G 

12  6 

2 28 

1000 

0 0312 

14  1 

3 35 

2000 

0 0180 

13  8 

3 03 

O,  (l 

r=  12  1 oV) 

60 

0 581 

119 

-20-7 

100 

0 411 

13  9 

-10  3 

200 

0 120 

16  7 

- 2-19 

300 

0 0S30 

10  9 

- 0-519 

600 

0 0093 

17  3 

- 0 295 

1000 

0 02S9 

1 9-2 

0 03S 

2000 

0 0190 

19  2 

• 0 4C0 

* A(E)  is  in  units  of  10~'e  cin*/oV,  \'{K)  nnd  T0(E) 
ore  in  oV. 

■J  Theso  values  nro  obtnined  by  fitting  Oi*al  el  al. 
(1971)  dntn  using  equation  (4)  with  k = 1. 

(Ref  3:1724) 


Table . A- 2 

Pnrnmeters  of  He,  Xs  and  O, 


r.(oV) 

r.(oV) 

A’,, 

A,(cV) 

K. 

J.(  oV) 

Ho 

-2-25 

16-5 

1-48 

0-34 

1-30 

1 -32 

N, 

4-71 

13  8 

7 12 

0 82 

0 85 

1 74 

o. 

1-86 

18-5 

7-90 

11-1 

7-07 

3-70 

• Tt  ~ 1000  cV,  Tb  >-»  21,  mid  1%  “ I.  Krv  nndJC|>  nroobtuinod 
by  filling  T)  of  Or.u.t!  al.  (1971),  and  K„  nnd  J,  nro  obtnined 
by  filling  cr ,(/>’)  of  and  llUN.v  (1000). 


(Ref  3:1726) 


71 


Loss  Functions  and  Reaction  Efficiencies 

Figures  A-4  through  A-7  illustrate  loss  functions  and  reaction  effi- 
ciencies obtained  by  Peterson  and  Allen  (Ref  6)  for  the  element  argon,  and 
Bass,  Berg,  and  Green  (Ref  4)  for  the  element  mercury. 

The  loss  function  is  defined  as  v 

S 

L(c)  = I to  .o  . (e) 
j 3 J 

+ E (Ij.0i (e>  + ;oE  1'>/  TSi(e,T)dT) 

where 

uk  = excitation  energy  of  the  jth  excited  state 
Oj  = excitation  cross  section  for  state  j 
1^  = ionization  potential  of  state  i 

= ionization  cross  section  for  state  i 
e = energy  of  the  incident  electron 
T = energy  of  ejected  secondary  electrons 

Reaction  efficiencies  are  defined  as  the  fraction  of  incident  energy 
going  into  each  of  the  inelastic  channels  and  can  be  expressed  as 

nj  = 

c 

where  J is  the  number  of  excitations  to  state  j produced  by  a primary 
electron  of  energy  e . 
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Mercury 


Fiyure  A-6 

Loss  (auction  and  cV/ion  pair  The  square  is  lire  eV/ion  pair  icsull  of  Jcs>:  (1971). 

(Ref  4:1863) 


I'flicicncie*.  iriine  Smith-  Lisla  ioniralion  cross  section 

(Ref  4:1863) 


W (£)  (cV/ien  peir) 


Appendix  B 


Justification  of  Initial  Assumptions 


Secondary  Electrons  Produced  via  a Mono-energetic 
Electron  Beam  Emerge  Isotrop ically  in  Velocity  Space 

Rigorous  analysis  of  the  velocity  space  isotropy  of  secondary  electrons 
produced  by  a mono-energetic  electron  beam  is  beyond  the  scope  of  this 
thesis.  However,  Mott  (Ref  2)  has  performed  cursory  examination  of  the 
phenomena.  Using  wavefunction  for  electrons  in  a given  energy  state,  moving 
in  a specified  direction,  in  a field  of  charge  Ze  ; and,  a wavefunction 
dependent  equation  for  differential  ionization  cross  sections,  he  obtained 
a relationship  between  the  energy  of  the  incident  electron  and  the  angle 
which  secondary  electrons  are  ejected.  Figure  B-l  illustrates  his  results. 


Angular  distribution*  of  scattered  and  ejected  electrons  corresponding  to 
electron*  of  200  volt*  incident  energy. 


(Ref  2:234) 
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For  incident  electrons  with  energy  much  greater  than  ejected  secondary 
electrons,  it  is  apparent  from  Figure  B-l  that  ejected  electrons  show  no 
real  directional  preference  and  hence,  can  be  considered  isotropic  in 
velocity  space. 


Elastic  Collisions  Between  Electrons 
and  Heavy  Particles  Can  Be  Ignored 

Consider  the  collision  in  one  dimension  of  an  electron  of  mass  m with 
a heavy  particle  of  mass  M (assumed  stationary) . Conservation  of  momentum 
and  energy  equations  require 

mu  = mu^  + MV  (B-l) 


and 


1 2 _ 1 2 1 .2 
^mu  = -^nu^  + ^MV 


(B-2) 


where 

m = mass  of  the  electron 
M = mass  of  the  heavy  particle 
u = initial  velocity  of  the  electron 
u^  = velocity  of  the  electron  following  the  collision 
V = velocity  of  the  heavy  particle  following  the  collision 
From  equation  B-2 

m(u2  -u2)  = MV2  (B-3) 

or 

m(u  + u ) (u  - u ) = MV2  (B-4) 

1 1 

But  from  the  conservation  of  momentum  equation  (B-l) 

(u  - u ) = ” v (B-5) 

1 m 
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Substituting  equation  B-5  into  equation  B-4  and  simplifying  terms  yields 


(u  + i^)  = V 


(B-6) 


Adding  equation  B-5  and  B-6  and  simplifying  yields 


M v 
m 


or 


V = 2 m u 
M 


(B-7) 


(B-8) 


This  last  result  can  be  used  to  determine  the  amount  of  energy  lost  by 
the  electron  during  the  collision  process  (or  the  amount  of  energy  gained 
by  the  heavy  particle)  as  a function  of  the  initial  electron  energy. 


1 2 


m / 1 


■^MV  = 4 - - mu 
M \ 2 


(B-9) 


Thus,  the  energy  transferred  is  proportional  to  the  ratio  of  the  masses  of 
the  particles  colliding.  Therefore,  one  can  conclude  that  the  energy  lost 
from  an  electron  when  undergoing  a collision  with  a heavy  particle  is 
on  the  order  of  10“3  of  the  electron  energy. 

Since  the  energy  lost  by  an  electron  in  an  inelastic  collision  is 
typically  an  eV  or  greater,  changes  in  the  electron  distribution  function 
due  to  elastic  collisions  between  electrons  and  heavy  particles  will  be 
small  compared  to  changes  resulting  from  inelastic  collisions.  This,  of 
course,  assumes  that  the  average  electron  energy  is  large  enough  for  a sig- 
nificant number  of  electrons  to  have  energy  greater  than  the  inelastic 
threshold  energy. 


The  Dominant  Electron  Loss  Mechanism  is  Due  to  Attachment 

In  many  of  today's  gas  laser  systems,  a significant  electron  loss 
mechanism  results  from  chemical  reactions  in  which  a molecular  species 
reacts  with  a free  electron  to  form  a neutral  atom  and  an  ion.  In  the 
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case  of  a mercury  chloride  laser  this  reaction  becomes  (Ref  10) 

Cl2  + e”  Cl  + Cl” 

If  the  plasma  analysis  is  limited  to  early  in  the  gas  breakdown  pro- 
cess, such  that  the  number  density  of  ionized  species  is  small  compared 
to  the  number  density  of  neutral  species,  recombination  of  electrons  and 
ions  can  be  ignored. 

If  one  also  assumes  infinite  geometry,  such  that  finite  boundaries 
are  eliminated,  electrons  physically  leaving  the  plasma  can  be  neglected. 
With  these  restrictions  on  the  solution,  the  dominant  electron  loss  me- 


chanism comes  from  dissociative  attachment. 


Appendix  C 


' 4SK 


Derivation  of  Analytic  Solution 


From  the  modified  Boltzmann  equation  derived  in  Chapter  II  and 
assuming: 

1)  Steady  state 

2)  No  superelastic  collisions 

3)  Attachment  only  at  zero  energy 

4)  Secondary  electrons  resulting  from  non-electron-beam  collisions 
emerge  with  zero  energy 

5)  Momentum  transfer  cross  sections  (Q)  are  independent  of  the  inci- 
dent electron  energy 

6)  Single  species  which  exist  only  in  the  ground  and  singly  ionized 
state 

the  Boltzmann  equation  becomes 


(E/N)z  d 


u df 
du 


+ c£(u)  + a (u  + I) (u  + I)f(u  + I) 


3Q  du 

- a . (u)  (u)  f (u)  + 6 (u)  /q  a,  (u)  (u)  f (u)du  - 6 (u)  41  a (u)  (u)  f (u)  =0 

1 1 3.  a 

(c-i) 


Since  the  significance  of  many  of  the  above  terms  varies  in  different 
regions  of  the  energy  axis,  solutions  to  the  Boltzmann  equation  above  and 
below  the  ionization  threshold  energy  (I)  were  considered  separately  (see 
Figure  C-l) . 
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Figure  C-l  Regions  of  Solution 


Constraints 

Initial  constraints  imposed  on  the  solution  include  the  following: 

1)  Piecewise  continuity 
fi(I)  = f2(I) 

2)  The  function  and  its  derivative  must  vanish  at  infinity 

lim  f (u)  = o lim  df (u)  _ 0 

u-h»  u-h»  du 

3)  The  function  must  be  finite  at  u = o 

4)  Electrons  must  be  conserved 

c&(u)  + $ (u)^o^  (u) (u) f (u)du  = 6 (u) ^a0a (u) (u) f (u) 


Region  1:  (u  >_  I) 

Allowing  for  attachment  only  at  zero  energy  and  assuming  non-electron 
beam  generated  secondary  electrons  emerge  with  zero  energy,  the  Boltzmann 
equation  for  u greater  than  the  ionization  threshold  energy  becomes 


1 


(E/N) 

3Q 


du 


u d_f 
du 


+ c£(u)  + oi(u  + I) (u  + I)f (u  + I) 


-oi  (u)  (u)  f (u)  = o 


(C-2) 
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Approximating  the  inelastic  cross  section  (Ref  10) 

- 


a.  (u)  =o(l  - i) 


then 


uoi(u)  = oq(u  - I) 

(u  + I)  o , (u  + I)  =0  U 
1 o 


(C-3) 

(C-4) 

(C-5) 


where  I is  the  excitation  threshold  energy  and  a is  a constant  dependent 

o 

2 

upon  the  species  and  reaction  considered  (approximately  10  cm  ).  Sub- 

“16 

stituting  C-4  and  C-5  into  equation  C-2  yields 
2 


(E/Nr  d_ 
3Q  du 


u df 
du 


+ c$  (u)  + o uf  (u  + I) 
o 


-o  uf  (u)  + a If  (u)  = o 
o o 


(C-6) 


or 


(E/Nr  d_ 

3Q  du 


u d£ 
du 


+ c£(u)  + 0Qu[f(u  + I) 


f (u) 


+0  If  (u)  = o 
o 


(C-7) 


if  f (u)  is  not  a rapidly  varying  function  for  u >_  I than 


[f  (u  + I)  - f (u)  ] 


df 

du 


Substituting  this  approximation  into  equation  C-7  yields 


(E/N)  ^ d 
3Q  du 


u df 
du 


c0(u)  + a ul  df  + o If  = o 
+ o—o 

du 


(C-8) 


(C-9) 


One  is  then  left  with  a second  order  differential  equation.  However,  if 


(E/Nr  d 


U djf 

du 


3Q  du 

can  be  considered  small  compared  to  the  remaining  terms  in  equation  C-9, 
an  iterative  technique  can  be  employed  to  obtain  an  approximate  solution 
in  which  the  second  order  term  is  ignored.  Using  this  approximate  solution, 
the  second  order  term  can  be  expressed  explicitly  and  reinserted  into 
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equation  C-9  to  obtain  increasingly  accurate  solutions.  This  assumption 
must  be  checked  for  self  consistency  once  a final  form  for  f(u)  is 
obtained. 

Using  this  approach  and  ignoring  the  second  order  term,  equation  C-9 
becomes 

o Iu  dfQ  + a If  = -c£(u)  (C-10) 

o o o 

du 


or 


d (ufG)  = -c£(u) 

du  oQI  (C-ID 


which  can  be  solved  to  get 

f = -c  _1  f0  (u)  du  + c/ 
° aQI  u u 


(C-12) 


Approximating  the  secondary  source  spectrum  (Reference  Appendix  D)  by 


S(u)  = 


So 

2u 


uoSo 


V 2u2 


° < ii  < uo 


uo  I u 


(C-13) 


where 

SQ  = /“g(u)du 

u = excitation  threshold  energy  I 
o 

and  substituting  C-13  into  equation  C-12  yields 
fo  = I / ?aPo  du  + £' 

V u 2u2  u (C-14) 


f = cUqSq  1 c’ 
O — — 

2oqI  u2  u 


(C-15) 
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Using  this  approximate  form  for  f (u)  , the  second  order  term  in  equation 


C-9  can  be  written  as 


(E/N)' 


or 


3Q  du 


(E/N) 


3 Q du 


u dfQ 

= (E/N)2  d 

u d /cuoso  1 \ 

du 

3Q  du 

du  20oI 

L \ 

- 

/ J 

U dfo 

= (E/N)2  d 

- cuoS0  1 

du 

3Q  du 

1 

Q 

O 

M 

1 

(C-16) 


(C-17) 


(Note:  The  c'  term  has  been  dropped  to  prevent  introducing  two  constants 

of  integration  when  essentially  a first  order  differential  equation  is 
being  solved.) 

Substituting  equation  C-17  into  equation  C-9  yields 


(E/N)^  d_ 
3Q  du 


cuoSo  1_ 
°oI  u: 


+ cB (u)  + 0 ul  dfl  + a If  = o 
o o 1 


du 


(C-18) 


which,  as  before,  can  be  rewritten  as 

d(ufl)  = - cS( u)  + (E/N) 2 cuoSo  d 

3Q  (%I)2  du 


du 


o I 
o 


(C-19) 


This  is  a simple  first  order  equation  which  can  be  solved  to  obtain  the 

form  of  f(u)  for  u ^ I . 

f 1 (u)  = - 1 fcB  (u)  + (E/N)2  cuoso  1_  + c" 

u °oI  32  (o0I)2  u3  u (c_20) 

Again,  using  the  approximated  form  for  the  secondary  source  spectrum 

(equation  C-13)  equation  C-20  becomes 


f^u)  = cuoso  1 + (E/N) 2 cuoso  1 + c 


2o  I 


3Q  (°0l)2  u3 


(C-21 ) 


To  satisfy  the  initial  condition  of  normalization 

co  L 

J0 u^f (u)du  = o 

it  is  clear  that  c"  must  equal  zero.  Thus,  the  form  of  f (u)  for  u >_  I 


is 
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(C- 22) 


f (u)  = cuoSQ  1 + (E/N) 2 cuOSQ  1 

2V  u2  32  (°oI)2  u3 

It  should  be  noted  that 

lim  f(u)  = o and  lim  df  (u)  = o 

uk“  U‘+*°  du 

as  required  by  the  initial  constraints  on  allowed  forms  for  f (u)  . 


Region  2:  ( o < u I) 

As  in  Region  1,  if  we  allow  for  attachment  only  at  zero  energy  and 
assume  that  non-electron-beam  generated  secondary  electrons  emerge  with 
zero  energy,  the  Boltzmann  equation  for  u greater  than  zero,  but  less 
than  the  ionization  threshold  energy  becomes 


(E/N)2  d_ 
3 Q du 


U df 
du 


+ c£(u)  + ck  (u  + I)  (u  + I)f(u  + I) 


+ a . (u)  (u)  f (u)  = o 

l 


(C-23) 


But  the  ionization  cross  section  for  u < I equals  zero.  Therefore, 
equation  C-23  becomes 


(E/N) ^ d 


U df 

du 


3Q  du 
From  equation  C-5 

(U  + Do.(U  + I)  = °oU 

and  from  Region  1 (equation  C-22) 
f(u  + I)  = cuoso 


+ c£(u)  + o. (u  = I) (u  + I)f (u  + I)  = o 
1 


(C-24) 


2°oI  (u  + I)2 


1 + (E/N) 2 cuoSo  1 

32  ‘V*  (u  . I)3 


Substituting  these  values  for  f(  u + I)  and  (u  + l)o^(u  + I) 
equation  C-24  and  solving  for  the  second  order  term  yields 


(E/N)2 

d_ 

u df 

= - c£(u) -aQu 

cuoso  1 

3Q 

du 

du 

_ ^O1  (u  + I)2 

+ (E/N)  2 cuoSo 1 

32  (0oI)2  (U  + I)3 


into 


(C-25) 
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or 


d_ 

du 


u 5?jL 

du 


3Qc  5(u)  - 3Q  cuoSo 


(E/N) 2 
— cuo^O 


(E/N)' 


21  (u  + I)' 


oQl2  (u  + I)3 

Direct  solution  by  separation  of  variables,  unfortunately,  results 
obnoxious  form  for  f (u)  which  appears  to  be  unbounded  at  u = o 
fore,  a modified  WKB  method  was  used  to  express 
= g(u) 


d_ 

du 


u df 
du 


as 


f (u)  = 2/g (u) du  - ug (u)  + const 
(Reference  Appendix  D) . 

Using  this  approximation,  equation  C-26  becomes 


f(u)  = 


3Qc 


(E/N) 

32 


2/5 (u) du  - u5(u) 


(E/N) ‘ 
- cuoSo 


°oI‘ 


CUqSq 

21 

2;  _ 


du  - 


(u  + I) 
du  - 


(u  + I) 


(u  + I) 


(u  + I) 


(C-26) 
in  an 
. There- 


(C-27) 


(C-28) 


+ c' 


(C-29) 


Substituting  the  approximation  for  the  secondary  source  spectrum  (equation 
C-13)  into  equation  C-29  and  performing  the  necessary  integrations  yields 
the  solution  to  the  Boltzmann  equation  for  o < u 1 
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To  complete  the  solution  for  f (u)  , c' ' ' must  be  determined  and  manipu- 
lations necessary  to  meet  the  initial  normalization  requirement  must  be 
performed. 


Apply  Boundary  Condition  to  Evaluate  c' ’ ’ 

c' ' ' is  obtained  by  insisting  that  f (u)  is  piecewise  continuous  for 
all  u . 


f (u)  = f2<u) 


(C-31 


From  equation  C-22 

fl(u)  = CUoSo  1_  + (E/N) 2 cu0S0  1_ 

2 ,23 

2oqI  U 3Q  (ocI)  u 

Evaluating  f-^lu)  at  u = I yields 


fl(l) 


cuS 
o o + 


2V 


(E/N)2  CUoSo 

3 Q 2 5 

o I 


(C-32 


From  equation  C-30 
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V 

i 


(E/N)  2ur 


3Q  Cuoso 


(E/N) 


2 21 


21n (u  + I)  + 21 


u 


+ cuoso 

«ol2 


(u  + 1)  (u  + I)' 

.,2 


(u  +1)  (u  + I)2  (u  + I)3 


+ c' 


Evaluating  f 2 (u)  at  u = I,  equation  C-30  becomes 


Substituting  C-32  and  C-33  into  equation  f^(I)  = fjd)  and  solving  for 
c"'  yields 


(C-34 


Normalization 

To  complete  the  solution  for  f (u)  , the  normalization  requirement  must 


be  fulfilled.  An  effort  was  made  to  utilize  the  inherent  normalization 


.HV.imnu..,.  ■ 


in  the  definition  of  f (u)  . In  doing  so,  it  was  necessary  to  obtain  the 
steady  state  electron  density,  n , which  appears  in  the  definition  of  c 
(equation  2-11).  By  employing  conservation  of  electrons,  one  obtains  the 
following  expression  which  can  be  solved  for  ng  . 

c/~  $ (u) du  + /"  oi (u) (u) f (u)du  = fa/“  oauf(u)du  (C-35) 


Unfortunately,  the  second  integral,  representing  electron  sources  resulting 

from  indirect  electron  beam  and  field  driven  ionization,  becomes  infinite 

when  evaluated  between  the  limits  of  zero  to  nfinity.  This  is  due  in  part 

to  the  approximation  made  for  the  ionization  cross  section.  Where  for  a 

true  cross  section 

lim  (u)  = o 
u 

the  approximation  becomes 

lim  o . (u)  = constant 
ux*> 


Thus,  in  introducing  the  various  approximations  required  to  obtain  f (u) 
its  inherent  normalization  has  been  lost. 

As  a consequence  of  this,  the  standard  normalization  procedure  was 
employed  in  which  a constant  defined  by 


/ u“f(u)du  = c 
o norm 


(C- 36) 


was  determined.  The  normalized  f (u)  is  thus  1/c  f (u)  . It  should  be 

norm 

noted  that  this  process,  though  straightforward,  is  a bit  tedious.  Evalu- 
ating equation  C-36  yields 

c « u**f  _ (u)  du  + /"u^f.  (u)du  (C-37) 

norm  o 2 I 

Substituting  the  forms  for  f^(u)  and  f 2 (u)  (equations  C-22  and  C-33) 

equation  C-37  becomes 
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BH 


— . . 


I 


' " | 


'norm 


3g 

CJ°fl 

u3/2, 

(E/N)2 

2uo 

3Q 

cuoSo 

2 11  > 

(E/N)2 

21 

o 

cuoSo 

2 fl 

U1'2 

3ol2 

o — 

(u 

+ I) 

,,,  .1  1/2. 
c / u du 

o 

cuoso  /j  u 3/,2du 

2/  u ' ln(u  + I)du  + 21/  u du  -/  u 


1/2 


5/2 


du 


(u  + I) 


(u  + I ) ' 


du  - I/1  u1/2  du  + f1 

o =-  o 


5/2 


du 


(u  + I)‘ 


(u  + I) 


2o  1 
o 


+ <E/N)2  CUoS_o  /"  u-5/2du 

,2  1 


3Q  (ool)' 


(C-38) 

Using  the  following  approximation  for  the  logrithmic  term  and  equation  C-40 
to  solve  the  more  difficult  integrals 


ln(u  + I)  = ln2I  + _u  - _1 
21  2 


(C-39) 


m 

'*n 

X 


m-1 


(2n-m-l) aX 


+ ; Xm~2  dX 

(2n-m-l)a  xn 


- - f dv 

(2n-m-l) a 1 x 


(C-40) 


where  X * ax^  + b\  + c (Ref  14)  one  obtains  the  following  expression  for 


'norm 


M ... 


(C- 


Equations  C-22  and  C-30  along  with  the  normalization  factor  C-41,  yields 


the  final  solution  for  f (u)  for  o < u < ® 


Appendix  D 


Analysis  of  Approximations  Used  in  Analytic 
Solution  Derivation 


Approximation  fo£: 


=*  g(u) 


For  a differential  equation  of  the  form 


or 


iff 

du* 


u d^f  + d£  =■  g (u) 
.2  du 


(D-l) 


it  should  be  possible  to  express  f (u)  as  a polynomial  series  near  u = o 

for  a non-singular  f (u) 

2 

f =*  a + a,u  + a u . . • 
o 1 * 


where 

uf ' ' -+  o as  u -*•  o 

and 

f'  a^  as  u -*■  o 

If  uf''  vanishes  at  zero,  f(u)  can  be  obtained  through  successive  approxi 
mation  (Ref  9).  Therefore,  using  the  iterative  technique 
df  : g (u) 

du  (D-2) 

and 

u d2f  - u d_Tdf  1 • U 

du  du  du  I du  (D-3) 
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Substituting  equations  D-2  and  D-3  into  equation  D-l  and  solving  for  f (u) 
yields 

f(u)  5 2/g(u)du  - ug(u)  + constant 

Note:  One  should  check  for  self-consistency  by  verifying  that 

u d2f  <<  df 
du  du 

over  the  range  of  u for  which  the  approximation  is  used. 


Appr oximation:  S ( u ) 

A reasonably  accurate  analytic  model  for  the  secondary  source  spectrum 
generated  by  a mono-energetic  electron  beam  has  been  developed  by  Green 
and  Sawada  (Ref  3)  and  is  defined  as 


f*(u) 


Nj 


eb  A(e) 
e 


L (u  - tq)2  + r2 


where 


A (E) 


T = T 
lO  1 S 


oQ  K In  E 
E J 

T. 


(D-4) 

(D-5) 

(D-6) 


E + r. 


(D-7) 


and 


N 

^eb 

e 

E 

Ts'rs,rb' 


10-16cm 

energy  of  secondary  electrons  being  produced 
number  density  of  the  species  being  ionized 
current  density  of  the  electron  beam 
electron  charge 

energy  of  the  incident  electrons  (electron  beam) 


K,  J = constants  dependent  on  the  species  being  ionized 


= 1000 


T =■  21 

b 

k - 1 

where  I is  the  ionization  threshold  energy. 

Unfortunately,  this  form  proves  cumbersome  when  doing  the  integrations 
necessary  to  obtain  a partial  solution  to  the  Boltzmann  equation.  There- 
fore, the  above  Lorentzian  form  for  S(u)  is  approximated  by 
( 


2(u) 


, o < u < u 


u < u 
o — 


(D-8) 

With  this  approximation,  the  volumetric  increase  in  electrons  per  unit  time 
produced  by  the  electron  beam  can  be  obtained  by  integrating  S(u)  over 
all  u . Let  the  volumetric  increase  per  unit  time  be  defined  as 

SQ  = /*  S(u)du 

= c ./JJ°  du  + c_/“  1 du 

1 o 2 Uq  — 

u2 


or 


So  = cluo  + c2 


(D-9) 


If  one  insists  that  the  approximated  form  for  the  secondary  source  spectrum 
be  piecewise  continuous 


lim  $ (u)  = lim  £(u) 
o 


u-*u“  u-^i+ 


then 


ClUo  = c_2  * C1  " c_2 

uo  .2 


(D-10) 
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Substituting  equation  D-10  into  D-9  and  solving  for  C2  yields 


c2  3 


(D-ll) 


Then  substituting  D-ll  into  equation  D-10  to  obtain  yields 


c ^ “ uo^o 

V"  (D‘12) 

Therefore,  the  approximation  for  the  secondary  electron  source  spectrum 
becomes 

f S„ 


$(u)  = < 


o 

2uo 

uoso 


V 2u 


o < u < u0 


uo  1 u 


(D-l 3) 


To  simplify  boundary  conditions,  let  u equal  the  ionization  threshold 

o 

energy  I . 

SQ  was  determined  rather  heuristically  by  analyzing  the  source  spec- 
trum resulting  from  the  Lorentzian  form  of  S(u)  developed  by  Green  and 

Sawada  and  choosing  S such  that 

o 

/^°  £L(u)du  ■ /^°  £ft(u)du  (D-14) 

where  S^tu)  is  the  Lorentzian  form  and  SA(u)  is  the  approximation. 


Table  D-l 


Parameter*  for  \nrious  pnscs* 


1 

T. 

r. 

r. 

^9 

J. 

N. 

210 

-6  49 

24  3 

21  6 

2 19 

13  8 

Art 

15  7 

0 87 

6 92 

- 7 85 

9 30 

3-75 

Krf 

14  0 

3 90 

7 95 

-13-5 

18-7 

8 43 

Xof 

12  1 

3 81 

7 93 

-115 

22  8 

5 95 

H*t 

15  4 

1-87 

7 07 

- 7-7 

2 07 

0 o:>89 

CH4t 

13  0 

3 45 

7 06 

-12  5 

18  7 

7-59 

H.O 

12  6 

1 28 

12  8 

12  6 

2 78 

0-0309 

CO 

14  0 

2 03 

13  3 

14  0 

5 39 

0 :>G8 

C,H,f 

11-0 

1 37 

9 28 

- 5 8 

119 

0-503 

NO 

95 

-4  30 

104 

95  - 

118 

0 ICS 

CO, 

13  8 

-2  46 

12  3 

13  8 

10  6 

0-530 

• Fitted  to  S(K,  T)  of  Or»X  cl  al.  nt  E — 500  cV  nnd  to  «,(£)  of  the  compilation  of  Kn'lirn 
and  Dunn  (lyCO). 

| Fils  improved  if  r,  »*  I. 

(Ref  3 s 1726) 
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Figure  D-l  illustrates  the  Lorentzian  source  spectrum  (equation  D-4) 

produced  using  parameters  obtained  from  Table  D-l  for  the  element  argon. 

Overlayed  on  Figure  D-l  is  the  approximation  to  the  source  spectrum  defined 

by  equation  D-13  with  uQ  equal  to  the  ionization  threshold  energy  and 

SQ  chosen  such  that  equation  D-14  is  satisfied. 

From  Figure  D-4,  S then  becomes 
o 

r t 


S = Au  N 3eb 
o o 


2.0 


-20 


(D— 1 5 ) 


where  A is  chosen  to  be  2.2  . 
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